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ie INTRODUCTION 


Prior to 1959 analysis and design of control systems 
was primarily performed by classical methods. These 
methods are adequate when there is a single variable 
parameter and a single loop system. 

In modern control engineering, where multiloop systems 
with characteristic equations with two or more variable 
parameters usually result, the above mentioned design 
methods become in general very complex. The designer has 
to complete several designs in a trial and error process 
and to choose the most suitable solution and to make 
usually some arbitrary choices based on his experience. 

Methods for studying the effect of variations in 
parameters of a multivariable, multiloop control system 
have been developed by Mitrovic [1] , Siljak [2] , Thaler 
[3] , and other POSE OLS. All these methods are based 
on the same concept and their dependence on a computer is 
characteristic. 

Mitrovic's [1] method which was introduced in 1958, 
allowed study of the effect of variations of two variable 
parameters on the location of the characteristic roots of 
a system. This method consists of specifying as variables 


the two lowest order coefficients of the characteristic 
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equation. Using then the characteristic equation as a 
mapping function, constant zeta, omega and sigma curves 
are transformed from the s-plane into the variable co- 
efficient plane. This method permits the designer to 
adjust these variable coefficients so that the roots of the 
characteristic equation may be set at desired locations. 

In 1964, Siljak [2] extended Mitrovic's [1] , co- 
efficient plane method into the Parameter Plane method in 
which the two variable parameters may appear linearly in 
any or all coefficients of the characteristic equation. 
The parameter plane method provides the designer with a 
procedure for factoring a characteristic equation and 
displaying the results in a parameter plane diagram, and 
allows him to obtain information about system stability 
and how it is affected with parameter variations. 

The parameter plane method was further extended in 
1965 [4] to a general case in which the two variable 
parameters may appear nonlinearly in the coefficients of 
the characteristic equation. Another extension of the 
Parameter Plane method to include more than two variable 
parameters has been defined as the Parameter Space method 
S| 

When in general only two parameters (or coefficients) 
are involved, for any complex value of s or a pair of (o,4,) 
values two simultaneous equations in two unknowns are ob- 
tained by setting the real and imaginary part of the 


characteristic equation equal to zero. 


iS 





In 1967, Bowie [6] while working with the character- 
istic equation of an inertially stabilized vehicle on the 
parameter plane, found that a complete set of roots for 
the characteristic equation could not be determined in 
the real root area by using the existing parameter plane 
theory. He further investigated this problem and found 
that specific pairs of 6,v, values, make the above men- 
tioned system of equations singular, which thus corresponds 
to a single straight line in the parameter plane, called 
a "singular line". 

The singular line added a new dimension to the para- 
meter plane because it provides the designer a way to hold 
a root fixed while varying others, by varying the two 
parameter values. Another important characteristic of 
singular lines is also the fact that they can be drawn in 
a parameter plane diagram by hand, without needing a 
computer, as in the case of drawing constantZ orw, curves. 

Although an analytic study of the two previous men- 
tioned simultaneous equations, which are obtained by 
setting the real and imaginary part of the characteristic 
equation equal to zero, makes obvious the existence of 
singular lines under certain conditions, they were not 
found earlier due to the standard computer programming 
practices and limitations. 

Singular line theory is based on the special case of 
linearly dependent parameter plane solution equations. 


From the standpoint of the design engineer it was 
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desirable that this theory be investigated in order to 
establish a standard general set of rules which would 
provide a method to obtain systems with singular lines. 

In this thesis this task is accomplished and possible 
applications of the singular line theory in control system 


design are investigated. 
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| Fe ee SINGULAR LINE THEORY 


A. SINGULAR CONDITIONS 


Consider the characteristic equation: 


a 
K 
F(s)= zi ee at (2-1) 
wherea, {k=0, 1, 2,...,n} , is a function of two variable 


parametersaandB®, 1.e€., 4,=g, (a,8), and s is a point on 
the s plane which can be expressed in rectangular coordi 
nates as: | 
s= - 5% +54, 1-0 (2-2) 

If for a given pair of ¢5,,, values {where |¢,| <1 and 
Was 70}, equation (2-1) has an infinite number of «and 8 real valued 
pairs which satisfy it, then for this equation singular conditions 
exist. In this case the point of the s plane which is defined by the 
pair Coo Yns 1s a singular point, and the characteristic equation is 
Singular. 

The locus of all points ( a, 8) on the Parameter Plane which 
satisfy equation (2-1) in the case of singular conditions is called 
the singular line. Each such line corresponds to a specific singular 


point or in other words to a specific pair of 2 Wns Values. 


S? 
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B. SINGULAR LINE ANALYSIS IN THE LINEAR CASE 


A simple case, which was analyzed and which appears in 
many practical problems is the linear case. In this case 
the two variable parameters and 8 appear in the coefficients 


of the characteristic polynomial in the general form: 
d= BhatCy Btd, (2=5) 


where bye ,fe and de are real constants. 


lies Analysis by Using Parameter Plane Equations 


Consider the characteristic equation (2-1), where 
a. is defined by equation (2-3). By introducing equation 


(2-2), into equation (2-1), and equating real and imaginary 


parts to zero, it can be shown [6] that: 


Boa+, B = - D, (2-4) 
B, a HG. Oye et ae D, 
where 
BES. Kb ow K : ay 
ie Gia Upe(s)) =Bp = 1) Boy SL ETO (C8) 
K=0 —) 


C=) (Yee, wo SU Ce) Cae ¥ 1 WS Ca ees 
0) 


K=O SH 
nh K a 
— ! = K ; 
Diz) CDSa oS uy 2) Dee) (-1)Kaw KU Ce) 
K=0 ua! 


Rs, 





The functions U,(o) and Up.1 0) are Chebyshev 
functions of the second kind for which Ui )=0, SECS} aut? 
and in general: 
UL (o)= 26U,_1(5) - Up, _5(5) (2-6) 


The system (2-4) can be written in the Matrix form 


as: 
By Cy O ; -D, 
Bo Cy |B “Do 
or 
Ax =B (227) 


When the solution of the system (2-7) is considered 
there are the following cases: 
a. No Solution 
When Rank (A) = 1 and Rank (A, B) = 2, which 


implies that: 


we, 
-— 
€2 


ee (2-8) 
2 2 Do 


ee, 


then the system has no solution. 
b. Unique Solution 
When Rank (A) = 2 and Rank (A, B) = 2, which 


implies that: 


ee, 


D 
a = (2-9) 
=e 


4 
2) ‘2 
br [ke 


» | 
“th. 


2 


then the system has a unique solution. 
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c. An Infinite Number of Solutions 
When Rank (A) = 1 and Rank (A, B) = 1, which 


implies that 


oe) 


di! Stet Sper (2-10) 
then the system has an infinite number of solutions. There- 
fore in this case, since there are an infinite number of 
pairs of a,8 values which satisfy the system (2-4) and in 
consequence the characteristic equation (2-1), it implies 


that for this equation singular conditions exist. 


2 Analysis by Using Complex Variable Theory 


The characteristic equation (2-1) is considered 


again. This equation can be rearranged in the form: 

Pgs eO it 5 (5)) Pot fe iGs)) = 70 (2-11) 
where f(s), f(s) and f.(s) are polynomials of s. Intro- 
ducing equation (2-2) into equation (2-11), these poly- 
nomials will be in general a complex number with some 
magnitude and an argument. Thus they can be written as: 

f,(s) = 1, e 

: jé 

f(s) r, e 2 


93 


£.(s) tz 


Introducing the above equations into equation 
(et) yields: 


ty eJ®la + r? eJ82 Bs Tz eJ®3 = 0 (2-12) 
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Using the Euler's formula eJ” = cose + jsine 


> 


equation (2-12) becomes; 
r,{cos6,+jsine, Jurr,{cose ,+jsin8,} btr={ cos 6,+jsine, } = 0 
or 


{r,cosd, %r,cos 6,8 +r,cosd , }+j{1r,sind jatr,sind,B+r sind, } = Q 


Equating real and imaginary parts to zero yields 
the system: 
r,Sin0,a+r,sin6, B=-r,sSinez 
(2-13) 
rj cos8,a+trcos@, B =-TzCOS6 z 
By applying the same reasoning as in paragraph one, 


the system (2-13) will have: 


a. No Solution 


If and only if: 


r,Sino, : r,S1n6 5 r,Sind ~ 
r,cos6, 1,C0s9, rCOsO 2 
or 


tan6, = tané, 7 tano. 
b. Unique Solution 


land “only 1£- 


tan6, # tané, # tan . 


c. An Infinite Number of Solutions 
If and only if: 


tano, = tano, = tano (2-14) 


Hpi 





Equation (2-14) implies that for the character- 
istic equation (2-1) which has been rearranged in the 


form (2-11) singular conditions exist, if and only if: 
oT = 6, + k,l 


0, = 8, + kot 


£,(s) a) 3(s) BD kil 
(2-15) 
f(s) =/f£2(s) + kot 


where k, and k, are integers. 


OF 


Se Singular Conditions for the Linear Case 


Summarizing previously derived results, the con- 
clusion is that, in the linear case where the character- 


istic equation is written in the general form: 


n 
F(s) -)° ay = 0 
k=0 


where a, =B, at Cc. B+ di; singular conditions exist if and 


OnLy ii: 


When the characteristic equation is considered 


in the general form: 


H 
| 


f;(s)a + f2(s)8 + £,(s) 


LD 





then singular conditions exist, if and only if: 


/£, (s) =/t3(s) + ki 
yl) VB ee + k,l 


4. Equation of a Singular Line 





When for a given pair off, Ws values singular 


conditions exist for the characteristic equation: 


Nn 
F(s) -)L ay. sk = 9 
k=0 


or 
£,(s)o + f(s) 8 a £.(s) = 0 


then, according to the previous analysis the two para- 
meters and§8 will be related by one of the following 


equations, respectively: 


B= - ?1 e PL Q= - va : <2 Ol (2-16) 
a ae 
or 
B = De see ig 2 piel « (2-17) 
r,sino. r5sim 5 r,Cos6 4 r,cosd, 


Each of these equations, is the equation of the 


) 


Singular line corresponding to the singular pair ca ire 
which was considered. This singular line is a straight 


line on the parameter plane which is tangent to the constant 
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tC, curve at the singular point LO; altace tound in this 
thesis, that this is also tangent to the constant Wi. curve 


at the same point. 


Di Methods to Test if a Point s is Singular 


When the characteristic equation is given and it 
is desired to determine if a given point (Zz, Wad» Cg a! 


and a0) 4 is singular the following methods can be used. 


a. Algebraic Methods 


By using the, W, values of the point under 
consideration, the coefficients of equations (2-5) can be 
evaluated. Then the given point will be singular if and 
only if equations (2-10) are satisfied. 

A second algebraic method consists of the 
evaluation of the arguments appearing in equation (2-15), 
after writing the characteristic equation in the form of 
equation (2-11). Then the given point is singular if and 


only if equations (2-15) are satisfied. 
b. Graphical Method 


By this method the characteristic equation 
must be written in the form of equation (2-11) and the roots 
of the polynomials f,(s), £,(s) and £,(s) must be evaluated. 
After these roots have been located on the s-plane, with 
the use of a spirule the arguments of the previous men- 


tioned polynomials for the point under consideration are 
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determined. Then this point will be singular if and only 


if, equations (2-15) are satisfied. 


c. Use of a Computer 


There are in general many ways by which a 
computer can be programmed to answer this problem. One 
such program is the <<Singular Point proeean >> which was 
developed with this thesis and which appears in Appendix 
H. In order for this program to be used in the case under 
consideration, the corresponding value of c of the point 
which is considered is used as an input data. Then at the 
computer print output the corresponding value of Ww, is 
checked if it is such that it makes the system (2-4) have 
an infinite number of solutions. When this is the case, 
the point ( 34, oe) which was considered is a singular point, 


otherwise it is not singular. 


6. Methods to Test if a Characteristic Equation is 
Singular | 


In the previous paragraph a characteristic equation 
was considered to be given and it was desired to determine 
if a point of the s-plane was a singular point for this 
equation. In this paragraph the problem is reversed, l.e., 
a characteristic equation is considered to be given and it 
is desired to determine if there is at least one point of 
the s-plane which is singular. For this problem the fol- 


lowing methods can be applied: 
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a. Analytic Method 


By this method the characteristic equation is 
written in the form of equation (2-11). Then, the argu- 


ments of f(s), f,(s) and £,(s) are evaluated as functions 


of © and Be ce — aoa 
f, (Cs) = 81( o,0_) 
£,(s) = 8,( o,0,) (2-18) 


f.(s) = 820 3,u,) 


Introducing equations (2-18) into the system 


(2-15) yields: 


816 Ww) z= g2(% Wd <3 kit 


(Zo) 
£,(,4,) = 82,0) + kot 


When the system (2-19) is solved for the un- 
knowns ¢ and Wo» any acceptable solution, i.e., any solution 
where the values of & and uj are real and also |Z] < 1 
and a 0, will define a pair of 4, on values, which cor- 
respond to a singular point of the characteristic equation. 
If at least one such point is found then the character- 


istic equation is singular. 


b. Graphical Method 


By this method the characteristic equation 


which is considered must be rearranged again in the form 


Ze. 





of equation (2-11). Then the roots of the polynomials 
£,(s), f,(s) and £,(s) are evaluated and they are located 
on the corresponding points of the s-plane. By the use of 
a spirule a search on the s-plane must be performed for 
points which satisfy equations (2-15). This search can be 
performed either by moving on constant ¢ lines or constant 
We cycles. This method will in general be useful when 

the polynomials of s which are involved in equation (2-11) 
are of low order and it is desired to find out if the given 


characteristic equation is singular for a restricted area 


of the s-plane which is of interest. 
c. Use of a Computer 


The Singular Point Program which was mentioned 
before and which was developed for this purpose corresponds 
to one among the many ways by which a computer can be pro- 


grammed to answer this problem. 
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1 I ee. SINGULAR LINE THEORY AND CONTROL SYSTEMS 


A. DEFINITIONS 

For a control system with two parameters qgand Bg 
appearing in its characteristic equation the following 
definitions come as a consequence of the singular line 


theory which was previously developed. 
1. Singular System 


A control system is singular, if and only if, 


its characteristic equation is singular. 
2. Singular Characteristic Equation 


The characteristic equation of a control system 
is singular, if and only if, the system of the two 
Simultaneous equations in two unknowns, aand 8, which is 
obtained by setting s = - cut jw 1-54 into the character- 
istic equation and requiring the summation of reals and the 
Summation of imaginaries to go to zero is singular, for at 


least one pair of ¢, w. values, where (eles i “and ee ee 


n 


tS, 





5m Singular Root 


A root of the characteristic equation of a 
control system is singular if it makes the characteristic 


equation singular. 
4. Singular Point 


Each singular root corresponds to a specific 
point of the s plane. Such a point is called a singular 


point. 
B. SINGULAR LINE THEORY IN CONTROL SYSTEM DESIGN 


In general the analysis and design of any control 
system is based on an accepted model for the system under 
consideration. In the case of a linear control system of 
any order the no zero, two complex conjugate poles system 
model is used as a standard one. The accuracy of this 
approximation is tightly related to the dominant mode 
concept. 

At this point the use of the singular line theory 
seems to be of utmost importance since, if a system is 
singular the complex pair which provides the singularity 
remains fixed and the locations of the other poles may be 


adjusted to some extent because an infinite number of a 


30 





and 8 pairs is available for their adjustment. [In addition 
some other performance criteria may be satisfied by an 
appropriate parameter variation. 

All these cases are analyzed and further investigated 
for other probable applications of the singular line theory 
in control system design. 


Ge SINGULAR LINE THEORY AND CONTROL SYSTEMS 


1. Basic Requirements 


From all that has been stated in section II, it 
seems that in order for the singular line theory to be 
applied on control systems design and compensation the 
following basic requirements must be satisfied, i.e., the 
system under consideration must be singular, a singular 
root of the system characteristic equation must meet 
problem specifications when it is considered the dominant 
root, this root must be made dominant by adjusting the 
parameter ® and B values, and finally the two parameters 
a and 8B must be adjustable. 

. Therefore, the initial problem in general will be to 
select the structure of the system or the required 
compensator for a given system, in order for the above 
baste requirements to be meet. 

In the following paragraph a preliminary analysis of 


the form of the characteristic equation of a control system 


31 





was performed, as the first step in developing the re- 
quired structure of the control system to obtain the 


desired singular conditions. 


2. The Characteristic Equation of a Control System 


in the Linear Case 


In general the characteristic equation of a con- 


trol system in the linear case will be in the form: 


k 
F(s) = )- a, Ss) = 0 (3-1) 


where a, = by atc} B +d) and by» Cy and dy are real 
constants. Equation (3-1) can also be rearranged in the 


following general form: 
£,(s)% +£(s)8 +£,(s) = 0 C52) 


where f(s), £,(s) and f.(s) are polynomials of s of any 
order. Assuming that £,(s) and f,(s) are not identical to 
zero, and that £2(s) = 0 or £2(s) # 0, yields the following 


two cases: 
a. £.(s) # 0 


By assuming that £.(s) *# 0, dividing equation 


(3-2) by £,(s) yields: 








f(s), £208) , 


>+] = Q (3-3) 
£z(s) £z(s) 


or 
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Se) eve) Gea i 34) 


where Fj (s) and F,(s) are either a ratio of two polynomials 
of s or a polynomial of s. 

When the characteristic equation of a control 
system is considered to be in the form of equation (3-4), 
then the necessary and sufficient conditions for this 
equation to be singular, according to equations (2-15), 
are: 


F,(s) = k,l 
(3-5) 
7 (Ss) ss ko1 


where ky and K, are integers. 


[It must be noticed that according to the process 
which was followed in the derivation of equations (2-15), 
they are valid if either f,(s), £,(s) and f.(s) are con- 
sidered polynomials of s or ratios of two polynomials of s. 
Based on this remark and the fact that the argument of 
unity is zero or an even multiple of 1, equations (3-5) 
were derived from (2-15). 

Considering now the characteristic equation of a 
control system in the form of equation (3-4) the following 


Cases may appear: 
(1) Case I: F,(s) # F,(s). In this case a 


point s of the s-plane is singular for the charactereristic 
equation (3-4), if and only if, it does not lie on the real 


axis and satisfies equations (3-5). 
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In order to determine if a given point in such a 
case is singular or in order to find the singular points of 
the characteristic equation (3-4), if any, one of the 
methods which were described in the previous section can 
be used. [In addition to these methods the following method 
which is related to the Root locus concept can be applied. 

The characteristic equation of the control system 
is: 

F,(s) a+F,(s) B+] = 0 (3-6) 
Suppose that the Root loci of the following equations are 
drawn on the same drawing: 


A Fj (s) = -] 


B F,(s) = -] 
C F,(s) = -] 
D F,(s) = -] 


where A and B are Root locus variable parameters which take 
values from zero to plus infinity and C and D are Root 
locus variable parameters which take values from zero to 


minus infinity. 


Based on the angle relation on which each Root 


locus is drawn, all the points of each of the above 
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Root loci will satisfy the following equations respec- 


tively: 


(2n + 1)n 


die 
/F2(3) = (2n + l)0 


/F,(S) = 2nil 
/F,(S) =  2nil 


where n is an integer. Therefore, the singular points of 
the characteristic equation (3-6), 1.e., the points of the 
s-plane which satisfy equations (3-5) simultaneously, will 
be the points where either the first or the third Root 
locus intersect with the second or the fourth Root locus. 
From the above described method it can be con- 
cluded that in the case under consideration either there 
are no singular points or there is a finite number of such 
points. Therefore, when in general F,(s) # F,(s) the 
characteristic equation of the control system will not 
always be singular and if it is singular the probability of 
finding a singular point which meets problem specifications 


as the operating point of the system will be very low. 


Cece stl eel Gn bs). Inlethis case 
the characteristic equation (3-4) of the control system 
becomes: 


F,(s) (atg) +1 = 0 C52) 
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Thus, the singular conditions given by (3-5) become: 


/F{(S) = kj0 (3-8) 


As it will be shown in the following case, 
when there is only one angle relation to be satisfied, 
like that which is given by equation (3-8), then there 
will be in general an infinite number of singular points. 
Considering equations (2-5) and (2-16), the equation of 
the singular line of any singular point will be in the 


general form: 


a+B = C 
where C is a real constant. This implies that although 
in a singular case there will be an infinite number of 
pairs of a and 8 values which satisfy equation (3-7), the 
roots of this equation will be fixed for any of these 
pairs. For this reason this case can be considered as a 
trivial one under singular considerations and not useful 


for applications. 
(3) Case III: F,(s) = F,' (s). When F,(s) = 


Cs), where n is an integer, then the characteristic 


equation (3-4) of the control system becomes: 


F,(s) at+F,” (s) B+1 = 0 Esee) 


and the singular conditions given by (3-5) become: 


/F,(s) = ky (3-10) 
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i.e., in this case a point s of the s-plane is singular, 
if and only if it does not lie on the real axis and satis- 
fies equation (3-10). 

In order to determine if a point of the 
s-plane is singular for the characteristic equation (3-9) 
or in order to find the singular points of this equation 
one of the methods which were described in the previous 
section can be used. [In addition to these, the Root locus 
method which was introduced in the first case can be used 
as follows: 

Suppose that the Root loci of the following 


equations are drawn: 


A F,(s) -] 


=i 


B FL(s) 


where A and B are Root locus variable parameters which 
take values from zero to plus infinity and from zero to 
minus infinity, respectively. 

Based on the angle relation on which each 
Root locus is drawn, all the points of each Root locus 


will satisfy the following angle relations respectively: 


/ Fi (Ss) 
Ss 


where k is an integer. Therefore, all the points of these 


2k ale} i 


Z2kIl 


two Root loci will also satisfy equation (3-10). After 
these remarks it can be concluded that in the case under 


consideration there will be, in general, an infinite 
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number of singular points. Also the probability of finding 
a singular point which meets problem specifications as the 
operating point of the system will be relatively high. 

In the simplest case where n = 2, i.e., 


F,(s) = F,“(s), the characteristic equation (3-9) becomes: 


F, (s) a+F, *(s) g+1 = 0 Ge 


and the necessary and sufficient conditions for this equa- 


tion to be singular will be given again by equation: 


/F, (s) = kil (ee2)) 


where k is an integer. 


(4) Case IV: F,(s) = FCs). In this case 
interchanging the two parameters o and 8 in the character- 


istic equation (3-4) yields the previous case. 


All the rest of the cases, which may appear can 
be considered to belong in Case I, and can be analyzed in 


a Similar way. 


b. f(s) = 0 


By assuming that f,(s) = 0, equation (3-2) 
becomes: 


£2) era (S) 6 = (3-13) 
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Since f,(s) and f(s) are not equal to zero, dividing by 
f(s) yields: 
f(s) 





+1=0 C3514) 


Or 


F,(s) A + 1 = Q (3-15) 


The singular conditions for this case as it can 


be derived from equations (3-5) are given by the equation: 


/F(s) = ki (5216) 


where k is a real integer or zero. It can be shown by the 

Same procedure as was done in Case III in paragraph 2a(3), . 
that in this case there is in general an infinite number of 
Singular points. Considering equation (2-5) and (2-16) the 
equation of the singular line of any singular point will be 


in the form: 
“7 = (3-17) 


where C is a real constant. This implies that although in 

a singular case there will be an infinite number of pairs of 
a and 8 values which satisfy equation (3-15), the roots of 
this equation will be fixed for any of these pairs. For 
this reason this case also can be considered as a trivial 
one under singular considerations and not useful for 


applications. 
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Si, Summary and Conclusions 


In the preliminary analysis of the character- 
istic equation of a control system, which was performed in 
the previous paragraph, it was found that there is a form 
of the characteristic equation, which in general has an 


infinite number of singular points. This is the following: 


F, (s) atF, (s) B+1 = 0 (3-18) 


or in the simplest case where n = 2, the above character- 


istic equation becomes: 
F,(s) atF,“(s) +1 = 0 (3-19) 


where F, (s) *# 0, and it is either a polynomial of s or a 
ratio of two polynomials of s. 

Therefore equation (3-18) corresponds to the most 
favorable case for the application of singular line theory 
in control system design and compensation. 

The next step is to find the structure of the 
system (or the required compensator for a given system), 
which has a characteristic equation corresponding to the 
form of equation (3-18) or in the simplest case to equation 


(3-19). This is done in the following section. 
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IV. ONE INPUT ONE OUTPUT LINEAR FEEDBACK SYSTEMS 


A. PRELIMINARY ANALYSIS 


In figure (4-1) the block diagram of a unity feed- 
back control system is shown. Suppose that it is desired 
to compensate the system by making use of singular line 
theory. Then according to the results of the previous 
section, the problem is to find the compensation structure 
or structures which provide a system with a characteristic 
equation in the general form: 

eeMECC ee ES (e)a= 0 ae 
1.€., a singular system. 

Considering a compensation structure then the adjustable 
parameters are in general all the compensator parameters, 
1.e., the poles, zeros and gain values, and the gain of the 
plant. Since the poles and zeros of a compensator can be 
either real or complex, a preliminary analysis is needed 
to choose the two parametersa and £8 from among the gain 
values. 

In the following sections which include a lomtiicd 
analysis of the singular systems it is shown how one of 
these parameters, i.e., the parameter &, was introduced in 
the transfer function of a singular cascade compensator in 
such a way that its value affects its poles.- (See Section 


VI C and equation (6-18)). 
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In the preliminary analysis of the problem of finding 
the compensation structure or structures which provide a 
system with characteristic equation in the form of equation 
(4-1), 1.e., a singular system, the five most commonly used 


compensation structures were analyzed as follows: 


il First Compensation Structure 

The first compensation structure which was con- 
sidered is that of figure (4-2). 

The characteristic equation of this system is: 

1 + G(s) + G(s) H(s) = 0 (4-2) 

Comparison of equations (4-1) and (4-2) reveals 
that in order to obtain a singular system the following 
relations must be satisfied. 


G(s) a Fg (s) 


H(s) = — Fg(s) lta 
1.e., the feedback compensator must have poles and zeros 
identical with those of the plant, and the two parameters 
a and 8 must be related with the gain values as 1s shown in 
equations (4-3). Introducing equations (4-3) into equation 
(4-2) yields: 

1 +aFg(s) +6 Fg*(s) = 0 (4-4) 


which is the characteristic equation of a singular system. 
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ine second Compensation Structure 


The second compensation structure which was 
considered is that of figure (4-3). 
The characteristic equation of this system is: 


1+G(s) G(s) = 0 (4-5) 


Comparison of equations (4-1) and (4-5) reveals 
that in this case there is no way by which a singular system 
can be obtained by defining the two parameters a and 8 
between the plant and cascade compensator gain values. 

In this case it was also noticed that when G(s) 


and G(s) are defined as: 


(s+z,) (s+z,) eae (st+z_) 
G(s) = kFy (s) = k ac > Seer as 
s (st+p,) (st+p,) ... (stp,) 
and 
N 
a Ss (stp,)...(Stp_)+8k(s+z,)...(S+z_) 
G(s) = a+8 kKEg(s) = eee ee ee 


s\(s+p,) (stp,)...(s+p,) 


the characteristic equation of the.system becomes: 
eel 
SS (s)+8k Fy (s) = 0 


Since the argument of the second and third term in 


the above equation are not affected by the appearance of the 
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k and Re factors in these terms respectively, the above 
equation is also singular. 

Based on the fact that the cascade compensator 
which was thus introduced in order for a singular system to 
be obtained has the same number of poles at the origin as 
the plant, this singular case (like the singular case, 
which was analyzed in Section V) has the basic disadvantage 
that the type number of the plant is not preserved, except 


in the case where the plant is type zero. 


a2. Wiha Compensation Structure 


The third compensation structure which was con- 
Sidered is shown in figure (4-4). 
| Suppose that in general: 
G(s) = KF, (s) 
G(s) = k FCs) 
H{(s) = kis F.,(s) 
The characteristic equation of this system will 


be: 
1+G(s) G(s) + G(s) G(s) H(s) = 0 (4-6) 


Comparison of equations (4-6) and (4-1) reveals 
that in order to obtain a singular system the following 


relations must be satisfied: 
G(s) G(s) = kk. Fy (s) aS) = OF (s) FCS) (4-7) 


and 
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H(s) = ky Fy(s) = = Fy(s) FL Cs) (4-8) 


Introducing equations (4-7) and (4-8) into equa- 


tion (4-6) yields: 


iT 
© 


ltaFg(s) F.(s) + Fg (s) FL (s) (4-9) 


which is the characteristic equation of a singular system. 
This case can be considered (under singular consid- 
erations) identical with the first one of figure (4-2), 


when G(s) GCs) is assumed to be one transfer function. 


4. Fourth Compensation Structure 


The fourth compensation structure which was con- 
sidered is that of figure (4-5). 

The characteristic equation of this system will be 
in general: 


1+G(s) H(s) + G(s) G(s) = 0 (4-10) 


Comparison of equations (4-1) and (4-10) reveals 
that for a given plant i.e., given G(s) = KF, Cs); there is 
an infinite number of ways by which G(s) and H(s) can be 
defined in order for a singular system to be obtained. 
Based on the fact that certain types of feedback compensa- 
tion are in common use, e.g., velocity and acceleration 
feedback, only the following cases were selected to be 
analyzed as those which have the highest probability of 


application in practice. 
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a. H(s) = a 
When H(s) iS a gain adjustment which is 
defined to be the parameter a, then in order for a sin- 
gular system to be obtained the following should be the 
transfer functions of figure (4-5). | 
G(s) = kFg (s) 
G(s) = 8Fg(s) (4-11) 
H(s) = a 
Introducing equations (4-11) into (4-10) 
yrelds: 
ltakFg(s) + 8kFy“(s) = 0 (4-12) 
Because the factor k does not affect any argu- 


ment relation, equation (4-12) will correspond to a sin- 


gular system. 
b, Hs) = as 

When H(s) is defined to be a tachometer feed- 
back then in order for a singular system to be obtained 
the following should be the transfer functions of the 
blocks in figure (4-5). 

G(s) & kk Cs) 

=a ase : 
G(s) ae iy (s) (4-13) 


H(s) = as 
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Introducing equations (4-13) into equation 
(4-10) yields; 
ltaks F.(s) + 8ks* Fo Cs) a0 (4-14) 


which is the characteristic equation of a singular system. 
It should be noted that a basic steady state 
performance requirement of the system shown in figure (4-5) 
is that the cascade compensator G.(s) must transmit signal 
at zero frequency, 1.e., its transfer function must not 
have any zero at the origin. But in this case it can hap- 
pen, if and only if, Fy (s) has at least two poles at the 
origin, which implies that the plant must be at least type 


two. 
Cee anes) a= asé 


When H(s) is defined to be an acceleration feed- 
back, then in order for a singular system to be obtained 
the following should be the transfer functions of the blocks 


shown in figure (4-5). 
G(s) = kF, (s) 
G(s) = gst Fy (s) | (4-15) 
H(s) = nee 


Introducing equation (4-15) into equation (4-10) 


yields: 


4 


Laelees Fo(s) + ks Fy“ (s) = 0 (4-16) 


which is the characteristic equation of a singular system. 


47 





It should be also noted in this case that in 
order for the cascade compensator GCs) to transmit signal 
at zero frequency its transfer function must not have any 
zero at the origin, i.e., FCs) must have at least four 
poles at the origin, which implies that the plant must be 
at least type four. Because such a plant usually cannot be 


found in practice this case will not be further considered. 


5. Fifth Compensation Structure 


The last compensation structure which was con- 
sidered is that of figure (4-6). 

By considering the product G(s) G.7(s) as one 
transfer function then under singular considerations this 


case turns to be identical with the previous one. 
B.. SUMMARY 


In this section a preliminary analysis was performed 
on the most commonly used compensation structures in order 
to find which of these can provide a singular system, and 
under what conditions, i.e., a system with characteristic 
equation in the form of equation (4-1). 

The results of this analysis have indicated that from 
a theoretical point of view there are many ways by which a 


Singular system can be obtained by properly selecting the 
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compensation structure and by properly defining the com- 
pensators and the two parameters a and 8B. 

The conclusion also from this analysis was that three 
singular cases must be considered for a further analysis, 
because they may be useful and applicable in control 
system design and compensation. These are the following: 
The first singular case corresponds to the compensation 
structure of figure (4-2), where the feedback compensator 
and the two parameters a and g are defined by equation 
(4-3). The second singular case corresponds to the com- 
pensation structure of figure (4-5), where the compensators 
and the two parameters a and g are defined by equation 
(4-11). The third singular case corresponds again to the 
compensation structure of figure (4-5), where the com- 
pensators and the two parameters a and 8 are defined by 
equation (4-13), where the plant is considered to be at 
least type two. 

In the following three sections which refer to these 
three singular cases respectively, a detailed analysis is 


performed. 
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vie FIRST SINGULAR CASE ANALYSIS 


In the previous section the << first singular case >> 
was defined as that which corresponds to the compensation 
structure of figure (4-2), where the associated transfer 


functions and the two parameters a and 8 are defined as 


follows: 
G(s) = aF 5 ($) = ita Ai anon (5-1) 
Ses (0) eee (s+p,) 
and 
H(s) = & Fo (si) =£ ee ee Sibi Mamta! (5-2) 
* s\(s+p,)(s+p.) ..- (s+p,) 


A. SYSTEM OUTPUT AND ERROR ANALYSIS 


The system's transfer function is: 


C(s) OL F, (s) (os 
R(s) 1 +aF, (s) *BF, (s) 





OT. 


| 
DRL 2 Ned. i a 2 
s+ +0.S s+ St. SZ. 
sof ]es*Pi0 | |i P| |! j)*8[ |(st2,) 
} =1 ek 


k=1 W=1 


Mm 9] 
N o {| (St+Z..) ca (s+p,) 





T(s) = s 


(5-4) 
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The output of the system is; 

CUS eS NGS) IG) CSS): 
and the error of the system is: 

E(s) = R({s) {1 - T{s)} (5-6) 

When the final value theorem was applied to equations 
(5-5) and (5-6) the steady state output and error were 
evaluated, assuming that all the poles of the transfer 
function T(s) were on the left half of the s-plane. 

The input signals considered are those, which are 
the most commonly used, 1.e., step, ramp and parabolic 


input. 


se Plant Type Zero 


In this case where the plant is type Zero, i.e., 
N = 0, the factor a which appears in equation (5-4) be- 
comes unity. 

a. Step Input 

When a step input R(t) = Au(t) or R(s) = A/s 
was considered then: 
m n 
C(t) = C(t) = lim ie) = ees 
Ss S 


eo 50 n Mm nh Mm 
Ita? aa | Z ee oie 2 


a ha 1 phen K jzi J 


pk 





and 


n m 
A al |p 2+ z ae 
io) e=wcee cr) =alemes——- toi (s)}- a a a 
ie cc S+0o 
| lp sal [2 | Tp tal ey 
= i k q=2 Jk= 1 k . 


b. Ramp Input 


When a ramp input R(t) = Btu(t) 


or R(s) = B / 3° was considered then: 
= —~ 41: B 2 
St) = C.,(t) = lim s —; TMs) = « 
‘Emales S*O S 
and 
a _ 4: B 
e(t) = e,.(t) = lim s —-+{1 - T(s)} = 
tre S>O S 


c. Parabolic Input 


When a parabolic input R(t) = ct-u(t) 


eek (S) = Dee was considered then: 
C(t) = C..(t) = lim s “{ T(s) = = 

Tra : S*O S 

and 

e(t) = e,.(t) = lim s — (Gea IUCS Merce 
Cay S>O S 


ote Plant Type One 


In this case where the plant was considered to 


be type one, i.e., N = 1, the steady error and output of 


the system were found as follows: 


a2 
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1Oa.ae< 


and 


e(t) 


ax 
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e(t) 
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and 


e(t) 
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a. 


oS 


oo 


SS 


Se 
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(t) 


(t) 


(t) 


C(t) 


e 


pes 


(t) 


Step Input 


= lim s A Ts) = 0 
S>O 2 


= im s < €1-T(s)} = A 
S7*oO = 


Ramp Input 


= lim s =. T(s) = Bi SS 


S70 S = 2 
eee 


BT 
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= lim s {1-T(s)} == 
SO 


Parabolic Input 


er ZC Eo 
= lim s 3 T(s) = 


s7O 


= lim s ££ {1-T(s)} = « 
$s7O 5 


Plant Type Two 


When a type two plant was considered, i.e., N = 2, 


the steady state error and output of the system were found 


as follows: 
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e(t) 


t?- =< 
and 


e(t) 


t> « 


c(t) 


t-a 
and 


e(t) 


t+ 


C(t) 


e-> cc 


and 


e(t) 


toa 


step Input 


= lim s *T(s) = 0 
sto. 


= lim s <°{1-T(s)} = A 
S+0 


Ramp Input 


= lim s =, T(s) = 0 
S>O S 


B 


= lim s y, {1-T(s) } = « 


iF 7) 


Parabolic Input 





= lim s 4s {1-T(s)} = « 
S>O S 


Plant Type Three or Over Three 


In all such cases C., (t) = 0, and the error be- 


comes finite or infinity depending on the input signal 


which is considered. 
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on Remarks 


From the above results which are summarized in 
Table I of Appendix E, it can be seen that in general the 
type of the plant changes due to the form of the feedback 
path which is introduced in order for a singular system to 
be obtained. Therefore, compensation of a given system by 
this method yields a system which in general, does not 
respond to a given input as might be expected. 

Although some of the above cases which were 
analyzed may have a specific application in control sys- 
tem design and compensation, only the case where both 
output and error have finite values at steady state seems 
to be of a more general interest. 

The example which follows is referred to that 


case. 


B. EXAMPLE (5-1) 


Assume that in figure (4-2) the transfer functions of 


the plant and the feedback compensator are as follows: 


ais) = ar (s) = g———-—> —____—_ (S27) 
& (s+1)(st+2) (s+4) 
and 
. 8 Pe Cpeee ee (5-8) 
H(s) fo} Fg (s) &® (st+1)(s+2) (st4) 


oS 





then the characteristic equation of the compensated sys- 


tem is: 
1+ oF, (s) + BF, “(s) = 0 (BES 


mueatter SubStituting Equations (5-7) and (5-8) into (5-9) 


yields: 


Peas C77als’ +(2i2 = 10e)s” 4(308 +35a%a) s- + 


+ (224 + 500 + 68) s + (64 + 240+98) = 0 (5-10) 


By using the method of paragraph 6c, section II, the 
Singular points of the system were found for several values 
of t. They are listed in Appendix A. One of these sin- 
gular points, i.e., the singular point (o. = 0.60, ee 
3.1104) was further considered. Then from equation (2-6) 
the Chebyshev functions of the second kind were evaluated 


and from equation (2-5) the coefficients By; C, and Dy 


were found to be: 


By = -5.5085 
Cy = 0.6746 
D, = 41.7741 


Then from equation (2-16) the equation of the correspon- 
ding singular line was found to be: 


8B = - 61.9254 + 7.8693a (5-11) 
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With a computer program then, the roots of the 
characteristic equation (5-10) were found for several 
different parameter pair values, which satisfy equation 
(S-11). They are listed in Appendix B. By examining 
these roots it can be seen that the singular root and a 
real root are fixed and the other three roots are varying 
for different pair values. 

In figure (5-1) the parameter plane diagram of this 
Singular system is shown for some arbitrarily selected 
values of ct and O 

In figure (5-2) the parameter plane diagram of the 
same singular system is shown for one constant g and one 
constant w, curve, which correspond to the singular values 
of the previously considered singular point, i.e., ena 
0.60 and ed 3.1104. On the same diagram the singular 
line of this singular point is also shown. It was noticed 
that these two curves and the singular line are tangent to 
each other on the same point of the parameter plane. 

For a further analysis of this system and of any 
other system which is referred to the << first singular 
case >>, the same procedure which is described in the next 
section and which is referred to the ‘second singular case’, 


can be followed. 


Se 





Vi. SECOND SINGULAR CASE ANALYSIS 


In Section IV, the "second singular case" was defined 
as that which corresponds to the compensation structure of 
figure (6-1) where the associated transfer functions and 


the two parameters a and 8 are defined as follows: 


(So Ae CS t25). weal StO) 
G(s) = kFg (s) = > (6-1) 
s (s+p)(stp,) ... (stp) 
(StZ (st 2). eee (SZ 4.) 
G(s) = 8F,(s) = 8 S-——*>____"__ (6-2) 
. s (s+p,)(s+p,) --- (stp) 
H(s) =a (623) 
ae SYSTEM OUTPUT AND ERROR ANALYSIS 
The transfer function of the system is: 
ccs) _ SS) GQCs) _ 


T(s) = SONS 
MGS Veen hea G(s)G_(s) 


Introducing equations (6-1) through (6-3)‘into equa- 
tion (6-4) yields: 


m 


Bk |_| (s+Z.)* ae 
ee 


2N n 2 \o m m 2 
s° | | (stp ,) +aks | [Cs+P,)] |] Gst2,)+8k] |(stZ,) 
k= k=1 j=1 j= 
(Os), 
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The output of the system is: 
C(s) = R(s) TCs) (6-6) 
and the error of the system is: 
E(s) = R(s) - C(s) = R(s) {1-T(s)} (6-7) 


When the final value theorem was applied to equations 
(6-6) and (6-7) the steady state output and error values 
were evaluated, assuming that all the poles of the trans- 
fer function (6-5) lie on the left half plane. This 
analysis was based again on the three usually considered 


input signals, i.e., step, ramp and parabolic input. 


1. Plant Type Zero 


In this case where the plant was type Zero, i.e., 
N = 0, the steady state error and output of the system 


were as follows: 
a. Step Input 


When a step input R(t) = Au(t) or R(s) = A/s 


was considered then: 


ABk r1,.2 


met) = C._(t) = lims= T(s) = a 


.- " [VPetoai Ing! Tay sand 


Keer J 


je 





and 


e(t) 


ton 


or R(s) 


C(t) 


C7 


or R(s) 


C(t) 


t>« 


and 


e(t) 


Cro 


C 


C 


22 


Ss 


, (t) 


B/s° 


(t) 


(t) 


(t) 


= lim s“{1- T(s)}= = = = 7 
S70 7 
Py. +ak P ae 
c= k k=1 Key J 
Ramp Input 


When a ramp input R(t) = Btu(t) 


was considered then: 


= lim s > T(s) =< 


57:0 


EN lez 


= lim s "2 {1-T(s)} =< 


S+O 
Parabolic Input 


When a parabolic input R(t) = ct“u(t) 


sme | 57 was considered then: 


= lim s —% T(s) = & 
S*O 5” 
Ze = 
Selim oa {1-T(s)} = « 
S+O S 
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tap Plant Type One 


In this case where the plant was considered to 
be type one, i.e., N = 1, the steady state error and out- 


put of the system were found to be as follows: 
a. Step Input 


When a step input was considered then: 


C(t) = C__(t) = lims £T(s) = A 
ae. Ss coe S 
and 
e(t) = e._(t) = lim s & {1-T(s)} = 0 
tx a S>O = 
b. Ramp Input 
When a ramp input was considered then: 
n my Se B _ 
elt) = C.,(t) = lim s =9 T(s) =« 
Cerace S>O S 
and n 
pa | 
Me = 6. (t) = lim s =, {1-T(s)}« — + 
(Cage 2 S79 S B a5 
je. 


c. Parabolic Input 


When a parabolic input was considered then: 
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Cét) = C__(t) = lim s -% T(s) = « 
tr ae S>O 32 
and 

_ ee ii 7 
e(t) = e,,(t) = lim Sz {1-T(s)} = « 
t>= S 


3. Plant Type Two 


In this case where the plant was considered to 
be type two, i.e., N = 2 the steady state error and output 


of the system were evaluated as follows: 


a. Step Input 


When a step input was considered then: 


c(t) = C__(t) = lims ST(s) =A 
trex s7*O 
and 
et) eee ee (t) = Tims — {1-T(s)} = 0 
tx S*O 
b. Ramp Input 
When a ramp input was considered then: 

= - 1; B = 
E(t) C_,(t) = lim s -, T(s) = 
t+ S>O S 
and 

= es nie B = 
e(t) = e., (t) > 16th Sy {1-T(s)} = 0 
toe S*O S 
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c. Parabolic Input 


When a parabolic input was considered then: 


: 2C 
C(t) = e__ (t) = lim s -3 T(s) = « 
tra == S>oO 5° 
and : 
; 2c Co oe be 
ec) - ©..(t) = Tim s ete a eel ee 
ha 


A Remarks 


As a first step toward a further study of the 
"second singular case", an error and output analysis was 
performed. From these Pe aits which were obtained and 
which are summarized in Table II of Appendix E, it was 
noticed that in contradiction with the "first singular 
case'"’ the type of the plant is preserved in this case, 
where the compensation structure and the associated com- 
pensators of figure (6-1) are introduced in order for a 


Singular system to be obtained. 


B. EXAMPLE (6-1) 


In this example a type one plant was considered, with 
transfer function. 


al (6-8) 


aes KEL is) = 100 —<—§ ——___ 
& s(st1) (st+5) 
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Suppose that this plant is compensated by the com- 
pensation structure shown in figure (6-1), then according 
to equations (6-2) and (6-3) the transfer functions of the 


compensators are: 


G(s) = 8 Fo(s) = 8 ——=~—— (6-9) 
s(st+1)(s+5) 
and 
H(s) = a (6-10) 
The characteristic equation of the compensated system 
1S: 
1 + G(s) H(s) + G(s) G(s) = 0 (6-11) 


Introducing equations (6-8), (6-9) and (6-10) into equation 


(6-11) yields: 


6 


Ss + 12s° 


+ 46s" +(604100a)s” + (25+600a)s- 
+500 as+1008 = 0 (6-12) 


Using the computer method of paragraph 6c, Section II, 
the singular points of the system were found for several 
values of <. They are listed in Appendix C. 

One of these singular points was further considered, 
1.e., the singular point (3 = 0.40, Wass 0.97096). Then, 
corresponding to the damping ratio value, Chebyshev 
functions of the second kind were evaluated from equation 


(2-6) and from equations (2-5) and (2-16) the equation of 
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the singular line which corresponds to the singular point 


which was considered, was evaluated to be: 
B = - 4.924970 + 0.25924 (6-13) 


Considering the last two terms of equation (6-12) a 
necessary condition for stability according to Routh's 
first criterion is that a>0 and @>0. Since a and 8 
parameters are related by equation (6-13), it implies that 


a necessary condition for stability is: 
0< a < 0.05265 (6-14) 


By arbitrarily selecting a = 0.05 from the range of 
values defined by (6-14), the corresponding value of 8 
was found from equation (6-13) to be, 8 = 0.013. Then the 
roots of the characteristic equation (6-12) were evaluated, 
for this pair of parameter values, by using a computer 


program with the ZPOLR subroutine. These roots are: 


55 ae 0.3884 + j0.8899 So = - 0.3884 - j0.8899 
7 a 9.2141 oye = moe Ue? 
S$. = - U2 9270 Sp = - 0.0701 


1.e., one of the roots is the singular root (ca 0.40,0 6= 


0.97096) or s = - 0.3884 + j0.8899. 
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In figure (6-2). the parameter plane diagram of this 
Singular system is shown for some arbitrarily selected 


values of Zt and wae 
Cc. HARDWARE CONSIDERATIONS 


In figure (6-3) the block diagram of a unity feed- 
back system (plant) is shown. The open loop transfer 


function of this system is: 


G(s) = kF,(s) = k ake as ree ETE? (6-15) 


aaa i a ae 
s (stp, )(stp,) ... (s+p,) 


Compensation of this system according to the "second 
Singular case'' scheme, yields the compensated singular 
System of figure (6-4). This compensation structure will 
be in general complex since it introduces a cascade and 
a feedback compensator. Therefore, a simplified version 
may also be desired. This modification is shown in figure 
(6-5), where the feedback compensator (gain adjustment) 
has been taken from the plant and has been introduced at 
the cascade compensator G(s). Also the Root locus gain 
8 from G(s) has been set as a gain adjustment in front of 
the minor loop. 

The characteristic equation of the systems in figure 


(6-5) and (6-6) are respectively: 


1 + ak F,(s) + 8k 20s) = 0 (GailGy 


F 
§& 
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and 
1 +a F,(s) + 8k F,“(s) = 0 (6-17) 


Comparison of these equations reveals that both are 
Singular and they have the same singular points, but they 
will be in general different, except when k = 1, in which 
case they become identical. By this modification the 
whole compensation structure has been introduced in cas- 
cade with the plant, as is shown in figure (6-5). 

An equivalent system of figure (6-5) is shown in 
figure (6-6) where the whole compensation structure of 
figure (6-5) is replaced by an equivalent cascade com- 
pensator. 

The following considerations can now be stated con- 
cerning the compensation structures of figures (6-4) and 
(6-5) or (6-6). First both are singular. Second with 
the compensation structure of figure (6-5) or (6-6), the 
compensation problem has been simplified and better 
organized since it will be required to build a filter from 
distinct components or a filter in an integrated form 
which will be introduced in cascade with the plant. This 
is a simpler and in general a more realistic approach to 
the problem. In the case also where K = 1, i.e., the gain 
of the plant is equal to unity then these structures are 
also equivalent. When K # 1 then, although they will have 


the same singular points (because the singular conditions 
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to be satisfied are identical, in both cases), the 
corresponding equations of the singular lines will be 
different. This means that for the same singular point, 
different points of the parameter plane (a, 8), will be 
associated and also the non-singular characteristic roots, 
will be in general different. 

The cascade compensator G.(s) which was introduced 
in figure (6-6) has the transfer function: 


B (st21).(st+22.). 20. (S+Z,) (S525) 


G (s) St Oh a a ae aa ae 
S (st+pi)(st+p2) ..- (stp) ta(stz;) «.. (s*z,) 


where the subscripts s c stands for << Singular compen- 
sator >>. Assuming that a pair of a and 8 parameter pair 
values which satisfies problem specifications has been 
selected, then this compensator and its transfer function 
can be precisely determined. 

This modified singular compensation structure of 
figure (6-5) or (6-6) was further analyzed as the most 
applicable due to its simplicity under hardware considera- 
tions, in comparison with that of figure (6-4). It must 
be also noticed that ina See application, where a 
gain feedback adjustment at the plant is possible or 
desired, e.g., when it is desired to eliminate the effect. 
of some nonlinearities that the plant may have, then the 
initial compensation structure of figure (6-4) may be 


considered for application. 
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D. STABILITY ANALYSIS 


lis Introduction 


One of the basic requirements of any control 
System is stability. The stability concept and related 
definitions can be found in reference [7]. The first 
question which arises when stability considerations are 
made for the singular structure of figure (6-6) is: What 
effect has the cascade singular compensator which has 
been introduced, on the stability of the initial uncon- 
pensated system? That is, in addition to the fact that 
this compensator makes the system singular, does it 
improve or impair the stability of the initial system? 

In either case the question which seeks an answer is: How 
1s it possible from an initially stable or unstable 

System to obtain a stable singular system? The interpre- 
tations of curves and singular lines on a parameter plane 
diagram and determination of stability is not in general 

a simple process. 

Stability analysis by using the Root locus 
concept seems to be relatively easy and can give a more 
complete picture of the behavior of the system. The 
Stability analysis as it can be performed by using the 
Root locus concept for the case of the compensated singular 


system is described below. 
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Za Stability Analysis by the Root Locus Concept 


By substituting Fy (s) from equation (6-15) into 
equation (6-17) and after some manipulations the character- 


istic equation of the compensated singular system becomes: 


a (s+z,)(stz,) ...(s#z,)s"(s+p; ((s+pg) ...(s#p.)+8k(s#21) v0. 


(stz_)°+ ac (s+p,) (stp) ..+(5*PQ) -4 (6-19) 


Assuming that the singular points of the singular 
compensated system have been found and that one of them 
has been selected as the operating point of the system, 
then the equation of the corresponding singular line can 
be evaluated from equations (2-5) and (2-16). This 


equation will be in general in the form: 
Siete INolaels: 


where A and B are real constants. By substituting the 
parameter 8 from this equation (6-19) and after some 
manipulations, the characteristic equation of the singular 


compensated system becomes: 


2 2 2 
Ak(s+z, ) (Gaz) ...(S#Z_) +s (stp, ) (stp, ) +. (S#p_) (tz, )6+z, ) ++. (S#z,) | 
oo SE ee eee 


2 p BED 2 2 2 
Bk(st+z,) (StZ.) ..(S#Z,) +S Nos+p, ) (stp, ) ...(s+p,) 


(6-20) 
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Since A, B and k are known constants the numerator 
and denominator of the above equation can be factored, 
which in general yields: 

CSsoAN ll S.sr4 eae: (s+2., ) 


(eS ees ee (6-21) 
(Seq, )) (sed inssipe cae ) 


where ae and P. may be real or complex. By treating the 
parameter a as the Root locus variable parameter the Root 
locus of the compensated system can be drawn, either by 
the classical graphical methods or by the use of a com- 
puter. In the last case where a computer is used, usually 
it is possible to proceed directly from equation (6-19) to 
draw the Root locus of the compensated singular system. 
Since it is known that as the Root locus variable 
parameter a of equation (6-21) varies from zero to infinity 
the roots of the characteristic equation move on the Root 
locus segments, starting from the poles and terminating 
to zeros or at infinity when zeros are not available, and 
that the singular root and maybe one or more other roots 
of the characteristic equation (6-21) are fixed, indepen- 
dent of the variations of the parameter a, it implies that 
some zeros and poles of equation (6-21) have to be 
identical. Therefore, the numerator and denominator of 
equation (6-21) must consist of a number of common factors 
equal to the number of the fixed roots that the charac- 


teristic equation of the singular compensated system has. 
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This number is at least two, corresponding to the 
system's complex singular root. This fact makes the pro- 
cess of drawing the Root locus of the compensated 
Singular system by the classical graphical methods, much 
Simpler than it initially looks. Thus the Root locus of 
the singular compensated system consists of two or more 
fixed points of the s-plane plus the classical Root 

locus segments. 

Considering now system stability all the related 
questions can be answered from the Root locus study. For 
example, if any Root locus segment or isolated point lies 
entirely on the right half of the s-plane, then there is 
no pair of a and 8 parameter values for which the compen- 
sated system can be stable. If on the other hand all the 
Root locus segments lie on the left half of the s-plane 
as well as all the isolated points then all pairs of a 
and 8 parameter values will yield a stable system. Finally, 
if some Root locus segments cross the imaginary axis in one 
Or more points, then the range of the parameter a values 
for which the compensated singular system is stable can be 
determined. The corresponding values of the parameter 8 
are then found from the equation of the singular line. 

The following example is an application of the above 
procedure and concepts which were introduced for the case 


of the singular compensated system. 


WZ. 





ai Example (6-2) 
In this example the uncompensated system (Plant) 
of figure (6-3), was considered to have transfer function: 


(ee (Gn (Gazz) 
S s(st+l) (st5) 


The characteristic equation of this system is: 
s”> + 654 + 55 + k = 0 (6-23) 


The Root locus of this system is shown in figure (6-7). 
The gain at the stability limit is k=30. When this system 
is compensated by a singular cascade compensator the block 
diagram of the compensated system will be that which is 
shown in enone (6-5) or the equivalent of figure (6-6). 
The singular cascade compensator which was thus introduced 
has a transfer function: 

Coy hl, (6-24) 

an s°+6s“+5s+o 
which was obtained from equation (6-18). 
The ehracceris cic equation of the compensated 


Singular system is: 


1+G(s) G_.(s) = 0 (6-25) 
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Introducing equations (6-22) and (6-24) into equation 
(6-25) yields: 


5°4125> 


+46s°+(60+a)S°+(25+6a)s+5astke = 0 (6-26) 
Two values of the gainkwere arbitrarily selected, 

i.e., k=20 andk=35, by which the uncompensated system 

becomes stable and unstable respectively. These cases 


were separately analyzed as follows: 


a. Uncompensated System Stable, k= 20 


For this value of k, the uncompensated 
system is stable and the characteristic equation (6-26) 


of the compensated singular system becomes: 
5°+125°+46s'+(60+0)S"+(25+6a)s-+5as+208 = 0 (6-27) 


In figure (6-8), the parameter plane diagram 
of the compensated singular system is shown for arbitrarily 
selected values off and w, . By choosing g, =0.40 the 
corresponding value of w,, was found in the computer by using 
the <<Singular Point Program>> as wyog= 0.97096. Therefore, 


a singular root of the compensated singular system is: 


2 
Sea > 5s Uy, + Ga,.v I=co = - 0.5884 + 3 0.8899 (6-28) 
From equations (2-5) and 2-16), the equation of the 
corresponding singular line was found to be: 
B = 0.24621la - 1.21242 : (6-29) 


1.e., if the pair of aand8 parameter values satisfies 
equation (6-29), then the characteristic equation (6-27) 


of the compensated singular system will have a fixed complex 
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root at the point s = - 0.3884 +j 0.8899. 
| Considering now the coefficients of the zero 
and first power of s terms of the characteristic equation 
(6-27), a necessary condition for stability is that both 
these are greater than zero, which implies that%? 0 and 
B>0. Based on this observation the roots of the character- 
istic equation (6-27) were found for a specific range of the 
parameter aandg values and with each parameter pair 
Satisfying equation (6-29), by using the subroutine ZPOLR 
in a computer program. These roots are listed in Appendix 
D. From this list several general remarks can be made 
about the range of the aand 8parameters in the stable 
region and the variation of the nonsingular roots in this 
region, i.e., the stable region in terms of the parameter 
a value is bounded between the values of 5 and 40 approx- 
imately, and that within this region the variations of 
the nonsingular roots is such that the dominance and 
sensitivity concepts are considerably affected, for dif- 
ferent parameter values. 

Substituting the value of 8 from equation (6-29) 
into equation (6-27) the characteristic equation of the 


compensated singular system becomes: 


6 a 
s°+125°+46s "+ (60+a)s°+(25+6a)s“+50s+(4.92420 -24.2485) = 0 (6-30) 


Treating then the parameteroa as the Root locus variable 
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parameter, equation (6-30) can be written as: 


; SORE Re. (7) ae Ge ai) 


594.196 24465 4+605°425s° -24.2485 


or by factoring numerator and denominator yields: 


(s+5.2230) (s“+0.7768s + 0.9428) _ 
2, a ES ee ee ee ae 
(s+5.2230) (s2+0 .7768S+0 9428) (st4.7195) (s+1.8458) (s-0.5653) 
(6-32) 

The numerator and denominator of the left part 
of equation (6-52) have common factors corresponding to 
the real root s = -5.2230 and the complex singular root 
s = -0.3884 + 30.8899. Therefore, the characteristic 
equation of the compensated singular system has a fixed 
real root and a fixed complex root at these points of 
the s-plane. 

The Root locus of this compensated singular 
system, which was obtained from equation (6-30). By the 
use of the ROOTLO subroutine in a computer program is 
Shown in figure (6-9). 

At this point it must be noticed that, the 
remarks stated in paragraph two about the fixed roots 
that a singular system has, the analytical results shown 
in Appendix D, and the Root locus of figure (6-9) are 
in agreement. It must also be noticed that the fixed 
roots of the characteristic equation (6-30) are shown by 


dark cross marks in a circle, in figure (6-9). 
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b. Uncompensated System Unstable, k = 35 


When the gain of the uncompensated system 
(Plant) has the value of 35, then the uncompensated system 
is unstable since a pair of complex roots of the character- 
istic equation (6-23) lie on the right half of the s-plane. 

In this case the characteristic equation (6-26) 
of the compensated singular system shown in figure (6-6) 


becomes: 


5°+125°+465 "+ (60+0)5°+(25+6a)5-+Sast+ 352 = 0 (6-33) 


In figure (6-10) the parameter plane diagram 
of this compensated singular system is shown. 

For the same value of Gas in the previous 
case, 1.€.,0. = 0.40, the corresponding value of Wis Was 
found by using the "Singular Point Program" as Ww F9.97096. 
Therefore, a singular root of the compensated singular 
system is s = -0.3884 u j3 0.8899. The fact that the singular 
root remains the same as in the previous case, although the 
characteristic equation has been changed comes as a conse- 
quence of the singular line theory which was developed in 
section II and III, 1.e., if the characteristic equation 
of the compensated singular system is considered to be in 
the form of equation (6-17), then in both cases Fg(s) is 
the same and only the value of khas been changed. Because 
the singular points according to the singular line theory 


are the points of the s-plane, which satisfy equation: 


/Fg (s) =n I 
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where n is a real integer or zero, it implies that since 
Fg(s) is the same in both cases the singular points will 
remain the same. 

From equations (2-5) and (2-16) the equation of 
the corresponding singular line was found to be: 

B = 0.140690 - 0.69281 (6-34) 

For the same values of the parameter %, as those 
which were used in the previous case, the corresponding 
values of the parameter Rwere found from equation (6-34). 
For these formed pairs the roots of the characteristic 
Bevation (6-33) were evaluated by using the subroutine 
ZPOLR in a computer program. They are also listed in 
Appendix D. By comparing the roots of the characteristic 
equation (6-27) and (6-23). which are listed in Appendix D, 
it can be seen that they are identical for the same value 
of the parameter q. The reason for this fact is that when 
the singular point s = -0.3884 + 30.8899 or (C = 0.40, 
aes = 0.97096), is considered, the equation of the 
corresponding singular line as a function of k, according 
to equation (2-5) and (2-16) is: 

B = = (4.92420 - 24.2485) (B38) 
When the value of the parameter 8 from equation (6-35) is 
substituted in the general form of the characteristic 
equation of the compensated singular system, i.e., into 


equation (6-17) yields: 


lta F,(s)+ (4.92420 -24.2485) Fy(s) = 0 
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or after substituting F,(S) from equation (6-22) and 


making the necessary manipulations yields: 


Sie des (G0anes (2586a)\s—-Sas+ 


(4.92420 - 24.2485) = 0 (6-36) 


i.e., the characteristic equation of the compensated 
Singular system is independent of the value of k when the 
same singular point is considered. Therefore, for the 

Same singular point and for the same value of the parameter 
a the compensated Singular system will have the same 
characteristic roots independently of the value of k. As 

a consequence of these remarks not only for k = 20 but 

for any value of k >0 the Root locus of the compensated 
Singular system will be identical with the Root locus 

shown in figure (6-9) when the singular point s = -0.3884 + 


j 0.8899 is considered. 


4. Summary and Remarks on the Stability of the 


Compensated Singular System 


When the plant to be compensated by a singular 
cascade compensator has no zeros in its transfer function 
then the parameter Band the Plant's gain k will appear 
only in the zero power of s term of the characteristic 
equation of the compensated singular system as can be 
seen from equation (6-19). In addition, only the product 
of Band k will appear in this term, plus a constant when 


the plant is type zero. When a specific singular operating 
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point of the system is considered Cao) then the para- 
meters aand will be related by the equation of the 
Singular line which corresponds to the above mentioned 


Singular point. According to equation (2-5) and (2-16), 


the equation of the singular line will be in the form: 


g = e Ce Bon (6-37) 


where A and B are constants. From this equation it will 
be observed on the parameter plane that in the case under 
consideration when the value of k varies it will change 
the slope of the singular line which passes from the 
selected singular point. From this equation also it can be 
seen that the product 8k can be expressed as a function of 
the parameteraalone. By substituting this product in the 
characteristic equation (6-19) of the compensated singular 
system the whole equation will be expressed as a function 
of the parameter a. Thus the roots of the characteristic 
equation will depend on the value of aonly. 

When the plant to be compensated by a singular 
cascade compensator has one or more zeros in its transfer 
function then the parameter 8 and the plant's gaink will 
appear in three or more power of s terms of the character- 


istic equation of the compensated singular system as can 


be seen from equation (6-19). Because they will always 
appear together in a product form the coefficient C, of 
equation (2-5) can be written in the form: 
n : , 
CST (le Ce anees U;_,(c) (6-38) 
aPTe 
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where Ch is the constant coefficient of the 8k product, 

in the jpower of s term of equation (6-19). When a specific 
operating singular point of the compensated singular system 
is considered, i.e., ( Ce)» then the parameters a and 

B will be related by the equation of the corresponding 
Singular line of the above mentioned singular point. Then 
according to equation (2-5) and (io). the equation of the 
Singular line will be: 

B = z (A + Bia) (6-39) 
where A’ and B! are constants. By continuing with the same 
reasoning as in the previous case, the same remarks were 
found to be applied here. 

Therefore it can be stated, that in general the 
stability of a compensated singular system when a specific 
Singular point is considered, is independent of the value 
of the plant gain k and depends on the value of the para- 
metera®alone or in other words the stability of the compen- 
sated singular system can be interpreted in terms of the 
parameteraalone. The stable region also (if any) on the 
parameter plane will be bounded between constantoa lines. 
Thus when the compensated singular system is stable for 
some range of the parametera values, it will also be 
stable for the same range of the parameter a values and 
for any value of the plant gain k. Therefore, when stability 
analysis is performed by the Root locus method where the 
parameter ais used as the Root locus variable parameter 


then the value of k can be set to any desired value as will 


be dictated from other design factors. 
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cE. HOW TO OBTAIN A SINGULAR ROOT AT A DESIRED POINT 
OF THE S-PLANE 


1. Introduction 
Compensation of a system by using a singular 

cascade compensator will be in general useful if the 
compensated system provides a singular root where the 
specifications of the problem dictate. Therefore, when 
designing by using the singular line theory the next 
problem which seeks an answer is: How a euler root can 
be obtained, where the system's specifications dictate? 

The general objective in all cases will be to 
make use of the advantage that, by varying the parameter 
values the singular root remains fixed and as other roots 
are varying the system's performance may be moved toward 


the system's specifications. 


ae Selection of the Singular Point 
Assume that a plant with transfer function: 


(Size (Sti) ween CS PZ) , ; 
G(s) # KF, (s) = ko (6-40) 
s (stp,)(st+p.) ... (s+p_) 
1s going to be compensated by a "singular cascade compensa- 


tor' whose transfer function according to equation (6-18) is 
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| (stz )(s*z)) are (s+z.) 





G(s) = 8 . 
2* s\(s+p_)(s*p,) -.-(s+p,)t a(stz )(s+z) --.(Stz,) 


(6-41) 
then the characteristic equation of the compensated 


Singular system will be: 


1 + oF, (s) + Bk F,“(s) = 0 (6-42) 


According to equation (3-10) the singular points 
of the system will be the points of the s-plane which 


satisfy the following equation: 


/Fg(s) = kyl (6-43) 


where k, is a real integer or zero. Therefore, the location 
of the singular points on the e emillaric is determined only 

by the zeros and the poles of the plant, i.e., the singular 
points are a function of PCs) only, and in consequence 

they can be affected only by changing Fi(s)- 

In the case where among the singular points there 
is such a point which satisfies problem specifications as the 
operating point of the compensated singular system, then 
there is no need to change FCs). Considering now the 
case where such a singular point is not available then Py(s) 
has to be changed in order to obtain a singular point where 
it is dictated by the problem's specifications. The fol- 


lowing examples were selected for illustration of the design 
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procedure in the above case. 
a. Example (6-3) 
Assume the plant to be compensated has the 
transfer function: 


G(s) = kF,(s) = RO cal (6-44) 
s (stl) (st5) 


and that problem specifications, i.e., damping ratio, 
natural undamped frequency, settling time and peak 
overshoot require a dominant singular root at the 


point (¢. = 0.4, w _ = 0.97096) or s_ = - 0.3884 + j0.8899. 


In figure (6-11) the poles of Fg(s) have been 
plotted on the s-plane as reference points, in order to 
measure the angles involved. The desired singular point 
1s also located on the s-plane. By the use of a spirule 


it was found that: 


fE ae ae ty + $5) = ~ 180° 


therefore the point S. 1s actually a singular point since 
for this point equation (6-43) is satisfied for k, = -l. 

This point also was found to be a singular 
point by the use of the "Singular point program" in 
example (6-1). 

It must be noticed that another way to check 
if the point ©. = 0.4 ,w . =0.97096) is singular for the 
compensated singular system, is to draw the Root locus of 
the uncompensated system. Then if the desired operating 


point lies on the above Root locus (according to the 
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analysis which was developed in paragraph 2a(3), Section 
III), it is singular. Because this method is more 

laborious its use in general is not justified, except in the 
case where the desired operating point is not singular and 
some trade offs have to be done in order to decide if Fg(s) 
must be changed (by the appropriate initial compensator) 

or another operating point which lies on the Root locus of 
the uncompensated system can be selected as the operating 


point of the compensated singular system. 
b. Example (6-4) 


Assume now that the plant to be compensated 
has transfer function: 


s+5 (6-45) 


G(s) = kF p(s) = k ee a’ 
s(st+10) (s+20) 


and that problem specifications dictate for an operating 


singular point (€. =0.6,w Lye eon S.= -0.64) 0.8. As in 


ns 
the previous example the poles and zeros of Fo(s) were plot- 
ted on the s plane as shown in figure (6-12). Then by the 


use of a spirule the argument of F,(s) was evaluated and it 


was found that: 


[F (5) Se eC ae melee See Cex127¢e10\5°45°)) a—. 247° 
(6-46) 


Therefore the desired point is not singular since it does 
not satisfy equation (6-43). This implies that FCs) must 


be changed somehow in order for the desired point to become 
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singular. From a theoretical point of view it can be 
achieved as follows: 

(1) Introducing an Additional Zero. Suppose 
that an additional zero is introduced at the point s=-1.82. 
Then the angle contribution of this zero will be +67° as 
it can be seen in figure (6-13a). Then, by defining the 
transfer function of the new obtained system after the 


above zero was introduced G(s), yields: 


G_(s) = kF_(s) = k{(s+1.82) Fg(s)} 


p(s) = - 180° 


which implies that the desired operating point ex = 0.6, 


and 


Wist 1) is singular for the system whose transfer function 
1S: 


G(s) ae (s+5) (st+1.82) 


s*(s+10) (s+20) 


(2) Introducing an Additional Pole. Suppose 
that instead of a zero an additional pole is introduced 
at the point s=-0.58, then the angle contribution of this 
pole will be -113° as it is shown in figure (6-13b). Then 
by defining the transfer function of the new system after 


the pole was introduced, G(s) yields: 


E = 1 
BA 8) = k Fy (s) = k Sane: Fe(s)} 


IR) =~360° 


and 
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which implies that the desired operating point (6,-0.6, 


Wiss 1) is singular for the system with transfer function: 
s+ 
G(s) = k SS 
P s“(st+0.58)(s+10) (s+20) 


(3) Introducing an Additional Pole and Zero. 


Suppose now that an additional pole and a zero are 
introduced at the transfer function of the initially 
given plant. This pair of points can be chosen in such 


a way that their net angle contribution will make 


the desired singular point (o.= 0.60,8..= 1) actually 


ns 
Singular. For example if the pole is chosen at the point 
Ss = - 16 and the zero at the point s = - 1.66 as it is shown 
in figure (6-13c), then by defining the transfer function of 
the new system after this pole and zero have been introduced, 
Gps) yields: 
st+1.66 
s+16 


[Fg (s) = [Fy (5) oer Oe eee GC 


which implies that the desired operating point is singular 


Grp(s) = k FL (s) = k { Fy (s)} 


and 


for the system with transfer function: 
(st+5) (st+1.66) 
G(s) = k 
Zp 


s?(st10) (s+16) (s+20 
3. Summary and Remarks 


Considering in general an uncompensated plant with 
open loop transfer function G(s) = k Fo{s), where Fo(s) 


1S a ratio of any two polynomials of s, then the argument 
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of Fy Cs) evaluated at any point of the s-plane will always 
lie between the limits: 

k,I < / Fg (5) < (k,+1) I 
where k, can be any real integer or zero. Therefore, the 
maximum angle correction which PCs) may require in order 
to satisfy equation (6-43) for a desired point of the 
s-plane will be less than 180° , which implies that a 
Single zero or a single pole or a zero and a pole or finally 
any number of zeros and poles can be properly selected and 
introduced to Fy (Cs) in order to make the argument of Fy (s) 
satisfying the angle requirement for singular conditions, 
for a specific singular point, which is desired to be the 
Operating point of the compensated singular system. 

These poles and zeros in general may be real or 
complex. In consequence it implies that when such a change 
of F,(s) is required an initial cascade compensator must be 
introduced, before a "singular cascade compensator" is 
introduced. 

Since a compensator with only zeros actually can 
not be built, the case of introducing one or more zeros 
alone has to be ignored as unrealistic. Therefore, among 
the rest of the possibilities available an optimal selection 
of the initial cascade compensator has to be done. This 
selection can be based on noise considerations, problem 
specifications and the "dominance of the singular root” 
concept which is described in the next part of this section, 


1.e., except the fact that the Root locus of the uncompensated 
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system will be forced to pass through the point where the 
desired singular root is located, the shape also of the 
resulting Root locus must be such that it promises some 
range of dominancy of the singular root. (See Example 
mo. O)). 

It must also be noticed that as it will be shown 
in the following analysis when an initial cascade compensa- 
tor is required, then this requirement will actually affect 
only the transfer function of the singular cascade compensa- 
tor and neither the complexity of the rest of the design 
problem nor the compensation structure (See Example (6-6) 


and figure (6-28c)). Under these considerations an initial 


cascade compensator with only zero's may also be considered. 


Ee DOMINANCE OF THE SINGULAR ROOT 


1. Dominant Root Concept 


The concept of the dominant root has wide applica- 
tions in analysis and design of linear systems. Considering 
a system with a characteristic equation of order n, then 
there will be n characteristic roots for this system. If 
one of these roots dominates the transient response, then 
this root will be called dominant. When the dominant root 
of a system is a complex root, s =—tw +jW, Wares , then the 
overall response of the system for a given input may be 
closely approximated from the corresponding standard second 
order system graphs (i.e., the second order system whose 


the characteristic equation is s/ 


+2cu_stw“=0 and has no 
Zeros). 
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2. Dominance of the Singular Root 


When a singular root has been determined from 
the problem specifications, then this singular root may 
be available from the given plant or not. [In the last 
case an initial cascade compensator will be required, in 
order for the desired singular root to become available. 
The design process of this initial cascade compensator 
was previously illustrated. 

In either case it will also be desired for the 
above considered singular root to be dominant for one or 
more (depends from the specific application) parameter pair 
values. The design process in all possible cases toward 
the dominancy of the singular root will be illustrated by 


the following examples: 
3. Example (6-5) 


In this example the Plant of example (6-2) was 


again considered with the Root locus gaink = 20, i.e.: 


1 
ENCES ESE Me A esis) (6-47) 


Since the value of the plant's gain does not 
affect the location of the singular roots on the s-plane, 
as was shown in example (6-2), the plant under consideration 
when it is compensated by a "singular cascade compensator" 
will have the same singular roots with these which are 
listed in Appendix C for the example (6-1) and (6-2). Some 


of these roots were considered and their dominancy was 


investigated. 
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In figure (6-14a) and (6-14b) the compensated 
Singular system is shown, where in figure (6-14a) the 
singular compensator is considered to be in integrated 
form, and in figure (6-14b) the same singular compensator 
is considered to consist of distinct components. 

The singular cascade compensator which thus has 
been introduced according to equation (6-18) has transfer 
Function: 

G..Cs) = Sesser (0-48) 
The characteristic equation of the compensated 


Singular system was found to be: 


Pomii2s +465 '+(60+a)s°+(25+6a)s-+5as +208 = 0 (6-49) 


a. Singular Root (5. = 0.40,0, = 0.97096) 


Ss 


When the singular root (¢,= 0.40,0 55 0.97096) 


S 
was considered, then the problem which was investigated was: 
Are there any pairs of parametera andg values for which this 
Singular root becomes dominant? 


From equations (2-5) and (2-16) the equation 


of the corresponding singular line was found to be: 


mee U.24635la - 1.21242 Die510) 


Introducing equation (6-50) into equation 
(6-49) the characteristic equation of the compensated 


Singular system in terms of the parameter a becomes: 
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eee Maes - (ana sis (25160) 5-505 


(GeaiD 
+(4.92420 - 24.2485) = 0 


Treating then the parameteraas the Root locus 
variable parameter, the Root locus of the compensated 
singular system for the singular point under consideration 
was drawn in the computer. It is shown in figure (6-9). 
The compensated singular system has three fixed roots, 1.e., 
the singular root s.=-0.3884 +3 0.8899 and a real root 
s=-5.2230, which are shown with dark circle marks at the 
corresponding locations of figure (6-9). The compensated 
Singular system also has three other roots whose location 
varies on the corresponding Root locus segments according 
to the value of the parameter a. 

Examining the roots of the characteristic equation 
(6-51) in relation with the Root locus of figure (6-9) 
it can be seen that the value of the plant gain does not 
affect the location of the roots of the characteristic 
equation of the compensated singular system and that the 
location of the three varying roots is affected only by 
the value of the parameter a. The dominancy also of the 
Singular root can be affected only by the other two 
Characteristic roots which lie on the same Root locus seg- 
ments with the singular root. These roots are either real 


Or complex according to the value of the parametera. 
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It was found that in the stable region the value 
of the parameter 4is approximately betwen the values of 
five and fourty. When a= 9.85 then the two varying roots 
in the vicinity of the singular root become identical with 
the singular root, 1.e., for the value of the parameter 
a = 9.85 the singular root is repeated. It was also found 
that a=9.85 not only the singular root but and the real 
fixed root at S=-5.2230 was also repeated. 

In an attempt to find the values of the parameter 
a for which the singular root is dominant (or in other 
words the system's performance can be well predicted from 
the corresponding to the singular root second order model) 
several runs with different parameter pair values each, 
were performed in the computer using a DSL (Digital 
Simulation Language) program. In this program the compensa- 
ted singular system shown in figure (6-14b) was implemented. 
When the second order model corresponding to the 
Singular root (Z. =0.40,w _ =0.97096 was considered with 


transfer function: 


f(s) _ cn : 0.9428 | (6-52) 
s+25.w s+ u s +0.77685s+0.9428 

its transient response characteristic values were calculated. 
These values together with those which were obtained from 
the computer for the compensated singular system are 


listed in Table I of Appendix F. 
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By examining this table the following remarks can be stated: 


(1) There is no pair of parameter values for 
which all the transient response characteristics of the 
compensated singular system can be well predicted from 
the values corresponding to the singular root second order 
model. 

(2) For all the parameter pairs both the 
Singular root and the other root in the vicinity of the 
Singular root contribute to the transient response of the 
compensated singular system. 

(3) Considering that in practice the only 
transient response characteristic values for which the 
designer 1s in general interested, are the settling time 
and the percent overshoot, then from this point of view, 
it can be said that the singular root is dominant for 
a =7.8 andg= 0.708. Assuming also that there is some 
tolerance in the system performance then these two parameter 
pair values can vary within some range according to the 
tolerance values. 

The conclusion from this example was 
that when dominant mode design is used for compensation of 
the given plant then the singular line theory provides in 
general an acceptable solution. In figure (6-15) the time 
response of the compensated singular system is shown for 
a unit step input and for the two pEemeon having the 
values of a= 7.8 and g= 0.71. In figure (6-16) the time 


response of the compensated singular system is shown for 
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the parameter pair values a= 9.85 andg = 1.21 for which the 
singular root and the fixed real root are repeated. This 

fact as can be seen from this figure causes a very high 
overshoot value. In figure (6-17) the time response of 

the compensated singular system is shown for such a pair 

of values, i.e.,a =20 and 8= 3.71, for which the contribution 
of the singular root and the other root in the vicinity of 

the singular root can be recognized. In figure (6-18) the 
time response of the compensated singular system is shown 

for the case where the characteristic root in the vicinity 

of the singular root lies on the real axis. In this case 
because this root is closer to the origin and to the imaginary 
axis than the singular root, it is dominant and the transient 
response of the compensated singular system has an exponential 


form. 


b. Singular Root G. = 0.30. =1.15268) 


The second singular root from the table of 
Appendix C which was examined for dominancy was the singular 
geoteG. = 0.3500, .= 1.15268) or s =-0.354580 = J 1.09959. 

From equations (2-5) and (2-16) the equation 
of the corresponding singular line was found to be: 
8B = 0.35266a - 2.48732 | (6-53) 

Introducing equation (6-53) into equation 


(6-49) the characteristic equation of the compensated 


Singular system in terms of the parameter abecomes: 


) 


6 ; 
$°+12s°+46s(60+a)5°+(25+6a)S-+5as+7.0532a - 49.7464 = 0 


(6-54) 
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Treating then the parametera as the Root locus 
variable parameter, the Root locus of the compensated 
singular system for the singular root under consideration 
was drawn in the computer. It is shown in figure (6-19). 
This Root locus as can be seen from this figure consists 
of three points of the s-plane where the three fixed roots 
of the compensated singular system lie and the classical 
Root locus segments. The above mentioned three fixed roots 
are shown with dark circle marks at the corresponding 


locations, 1.ée., 


S$, = - 5.3084 
Sy = - 0.34580+ j1.09959 
Sa, = 0.34580 -j51.09959 


In an attempt to find the values of the parameter 
a for which the singular root is dominant several runs with 
different parameter pair values each, were performed in the 
computer using the DSL simulation program. The transient 
response characteristic values which were obtained from 
the computer, together with these corresponding to the 
Singular root (ae =0.30,w, = 1.15268) second order model 
are listed in Table II of Appendix F. 

By examining this table, the same general remarks 
and conclusions were derived in this case as in the previous 
one. The basic remark again was that there is a pair of 
parameter values, i1.e., @= 10.6, B= 1.25, for which the 
Singular root is dominant from a practical point of view. 


In figure (6-20) the time response of the compensated 
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singular system is shown for this pair of parameter values. 
In figure also (6-21) through (6-23), the time response of 
the compensated singular system is shown for several 


characteristic pairs of parameter values. 


c. Singular Root (T = 0.75,W, 55 0.60667) 
The third singular root which was examined 


for dominancy was ( . = 0.75, 0.60667) or 


ns 
s = -0.45500 + 30.40127. 
From equations (2-5) and (2-16) the equation 
of the corresponding singular line was found to be: 
B = 0.093674 - 0.17547 (6-55) 
Introducing equation (6-55) into equation 


(6-49) the characteristic equation of the compensated 


Singular system in terms of the parameter a becomes: 


Mestas +465 +(60+a)s "+ (25+60) s-+5as+1.87340-3.5094 = 0 (6-56) 


Treating then the parameter aas the Root locus 
variable parameter, the Root locus of the compensated 
singular system for the singular root under consideration 
was drawn in the computer. It is shown in figure (6-24). 
This Root locus (as can be seen from figure (6-24) also 
fensists of three points of the s-plane where the three 
fixed roots of the compensated singular system lie and the 
Classical Root locus segments. The three points where the 
three fixed roots are located are shown with dark circle 


marks at the corresponding locations, 1.e.: 
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Sy= - 5.0899 


S55 - 0.45500+ j0.40127 
Sz— - 0.45500- j0.40127 

By the same procedure used in the previous cases 

the same general remarks and conclusions were derived in 
this case also. The basic remark again was that there is 

a pair of parameter values aand 8 for which the singular 

root is dominant under a practical point of view. This pair 
as can be seen from the Table III of Appendix F is: = 3.6, 
8B = 0.162. In this table the transient response character- 
istic values which were obtained from the computer together 
with these corresponding to the singular root (o. =tie! oO 
Whe = 0.60667) second order model are listed. In figure 
(6-25) the time response of the compensated singular sys- 
tem is shown for the previous mentioned parameter pair 


values. 


d. Remarks 


It must be noted here that the proceeding 
design was available for the particular system. It may 
not be generally available or may not provide an accept- 


able solution to some problems. 


4. Root Loci Correlation 


In the previous example if the desired singular 


root of the compensated singular system could not be made 
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dominant by a pair of parameter values, then an initial 
cascade compensator would be required, as has been already 
stated in paragraph two. 

Since the initial cascade compensator is applied 
to the uncompensated system, i.e., the system before the 
Singular cascade compensator has been introduced, in order 
to affect the compensated singular system, i.e., the 
system after the singular cascade compensator has been 
introduced, a prior knowledge of this effect is necessary 
in order to avoid an undetermined number of trials. 

This prior required knowledge or the criterion 
for the selection of the initial cascade compensator can 
be based on the correlation between the Root locus of the 
uncompensated system and the Root locus of the compensated 
Singular system, as can be seen from the following analysis. 

When an uncompensated linear control system with 
unity feedback is considered, its characteristic equation 
will in general have the form: 

(stz J(st253 ces (s+z_) . 


. , = k Fo(s) = - 1 (6-57) 
Ss (stp,)(stp,) ... (s+p,) 


Where k Fg(s) 1s the open loop transfer function and k 
is the Root locus gain. The Root locus of this uncompen- 
sated system which is drawn with the Root locus gain kK as 


the Root locus variable parameter consists of all the 


points of the s-plane which satisfy the angle requirement: 
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/ (s) = (2n+1)t 


Therefore any point Spr of the Root locus of the unconm- 
pensated system as it can be seen from equation (6-57) 
and the previous angle relation satisfies the following 


two relations: 


k PACs == 1 


RL? 


and 


/F g (Spy) = (2n+1)q 


or Fg(Sp1) < 0 (6-58) 


and /¥ g(Spy) = (2ntl) 1 (6-59) 


where n is a real integer. 


When the system which was considered is com- 
pensated by a singular cascade compensator the charac- 
teristic equation of the compensated singular system 


according to equation (6-17) will be: 


lta Fy(s) +8 k Fy “(s) =m (6-60) 


Then according to equation (3-12) a point of 


S-plane is a singular point for the compensated singular 
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system if and only if it does not lie on the real axis 


and satisfies equation: 


eo e'*) = nk (6-61) 


where n is a real integer or zero. Therefore all the 
points of the Root locus of the uncompensated system, which 
do not lie on the real axis (this restriction is implied 
from the definition of the singular point in Section II, 
A) are singular points because they satisfy equation (6-59) 
and in consequence equation (6-61). 

Suppose that one of these points, i.e., the 
point (¢.,w,.) is selected as the operating point for the 


compensated singular system. For this point equation 


(6-58) becomes: 
Fg(t.o4,,) <0 R626) 


The equation of the corresponding singular line 
(i.e., the equation which defines the relation between the 
two parameters a and 8 for which the selected singular 
point will always be a root of the characteristic equation) 
can be found either by using equations (2-5) and (2-16) 
or by direct substitution of the selected singular point 


in equation (6-60), which yields: 


Z : _ 2 
lt+a ele) eens Ce. 05) = 0 (66S) 
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Solving this equation for gk and substituting 


it in equation (6-60) yields: 


l+a F, (0 jw). 
$s" asp *(s) = 0 (6-64) 
2 


Fo (6.5, .) 


|e Fy cs) 7 


Equation (6-64) is the characteristic equation 
of the compensated singular system as a function of the 
parameter a, when the singular point (2.20, .) 1s considered 
to be the operating point of the compensated singular sys- 


tem. This equation can be written in the form: 


, ‘Fe 74s) Fg (S07 {Fg sens! ; Ba teyh © a. (6-65) 
{Fy (6554, 5) *Fy (s)} {Fy (E59 ¥, 5) - Pfs) } 


or 
(6-66) 


N 
Fy (6, 24,,)5 (stp, ) (stp)... (stp,)+(stz,) (stz)... (sz) 
Considering the Root locus of the uncompensated 
System, the Root locus of the compensated singular system 
and the corresponding equations, the following remarks 


and conclusions can be stated. 
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a. Number of Fixed Roots 
From equation (6-65), the common factor in 
the numerator and the denominator on the left side of 
equation, implies that the compensated singular system 
has a number of fixed roots equal to the order of the 
characteristic equation of the uncompensated system, e.g. 
when the order of the characteristic equation of the 
uncompensated system is three, there will be three fixed 
roots, two of which have to be the complex pair of the 
Singular roots. 
b. Root Loci Correlation 
The angle relation which the points of the 


Root locus of the compensated singular system satisfy is: 


(2n+1)1 (0507) 





Then according to equations (6-58) and (6-62) all the 
points of the Root locus of the uncompensated system satis- 
fy equation (6-67). Therefore all the points of the Root 
locus of the uncompensated system will be also points of 
the Root locus of the compensated singular system. 

Second, from equation (6-57) and (6-66) it 
1s observed that they have, at the expression on the left 
Side, identical zeros. They also have the same number of 


poles. 
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The conclusion from the above remarks is 
that the Root locus of the compensated singular system will 
be identical with the Root locus of the uncompensated sys- 
tem with the only difference that the Root locus segments 
of the compensated singular system may be further extended 
from the pole side end points, where the corresponding 
Root locus segments of the uncompensated system terminate, 
until they reach the poles of the compensated singular 
Syoctem. 

Examining the figures (6-7), (6-9), (6-21) 
and (6-25) they confirm all the above derived remarks and 
conclusions. Therefore, if an initial cascade compensator 
1S introduced to the uncompensated system this will affect 
the Root locus of the uncompensated system and the Root 
locus of the compensated singular system in the same way, 
i.e., the prior knowledge of the effect of the initial 
compensator on the compensated singular system can be 
obtained by studying the effect of the initial compensator 


on the uncompensated system. 


59. Range of Dominancy of the Singular Root 


When a singular compensated system is considered 
and an operating singular point, then the parameters a and 
8 will be related by the equation of the corresponding 
Singular line. In general there will be two ranges of 


a and 8 parameter values of primary interest. The first 
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one will be that for which the compensated singular system 
is stable and the second one will be that within the pre- 
vious one, in which the singular root is dominant. In 
general the wider these ranges the more desirable will be 
the solution. Therefore when an initial cascade compensator 
will be introduced to the uncompensated system in order to 
make the singular root dominant, it must be selected in 
such a way that it will also give a wide range of dominancy 
of the singular root. 

It must be realized that in a self-adaptive 
Singular system where the two parameters a and 8 are 
automatically adjusted so that the singular point will 
remain on the singular line during the operation of the 
system, the range of dominancy of the singular root will 
be very important. 

The following example is mainly referred to the 


range of dominancy of the singular root. 


6. Example (6-6) 


In this example the same plant which is shown 
in figure (6-26) was again considered. Its transfer 


function iS: 


age 
CD ie s(st+1) (st5) 
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For this example, it was assumed that problem specifica- 
tions require an operating point for the system, l.e., an 
operating dominant singular root, with damping ratio 

c = 0.60, and natural frequency Ww, = 1.2 rad/sec. A wide 
range of dominancy was also assumed to be required for the 
above singular root. 

The Root locus of the plant with the Root locus 
gain used as the Root locus variable parameter is shown in 
figure (6-29), where the location of the desired singular 
root is also shown. By studying this figure and taking 
into account the remarks which were stated in part VI E 
and VI F, the conclusion was that an initial cascade com- 
pensator is needed for the following-two reasons: The 
first is that since the desired singular root does not 
lie on the Root locus of the plant this Root locus must be 
reshaped in order to pass through the point where the 
desired singular root is located, and the second reason is 
that the shape of the Root locus must be changed in order 
for this singular root to become dominant with a wide range 
of dominancy, as it is required from the problem Specifi- 
cations. 

After the above remarks and taking into account 
the noise problem it was decided to compensate the plant 


initially by a cascade lead filter with transfer function: 


Si 
G(s) = S570 (6-68) 
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The compensated system is shown in figure (6-27). 


The characteristic equation of this system is: 


5 +4+265°+1255-+(100+k) s+2k = 0 (6-69) 


where k = 20. Equation (6-69) can be written in the form: 


= - id (6-70) 
s(s+1) (s+5) (s+20) 


The Root locus of this system with k used as the 
Root locus variable parameter is shown in figure (6-30). 

At z = 0.6 the corresponding Root locus point is 
me = 0.60, w= 1.232586). The value of W although it 
could be measured directly from the Root locus diagram of 
figure (6-30) was found in the indicated accuracy by the 
use of the "singular point program". Therefore this 
initially compensated plant provides the desired singular 
root (assuming that the small difference between the desired 
value of Wo» a= eae OL ES 1.2 rad/sec and the actual value 
of Wy = 1.232586 is within the system's tolerance). The 
Shape also of the Root locus of figure (6-30) is such that 
promises some range of dominancy of the singular root. 

Based on these remarks the initially compensated 
plant was next compensated by a singular cascade compen- 
sator the transfer function of which according to equation 
(6-18) is: 
ee 5) | (6-71) 


G 
se'8) ~ S(se1) (s+5) (5420) +a (5*2) 
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The compensated singular system is shown in 
figure (6-28a). The characteristic equation of this 


system was found to be; 


5°4525/+926s°+ (670040) 5" +(20825+28a) $*+ 


(25000+177a)s°+(10000+350a +208)s¢+ (6-72) 
(200a+808)s+808 = 0 


When the singular root (S. = 0.60, te 1.232586) 


was considered, the equation of the corresponding singular 


line was found from equations (2-5) and (2-16) to be: 
B = 3.31245248a - 219.4468287 (Gia -5)) 


Introducing equation (6-73) into equation (6-72) 


and after some manipulations yields: 


F 554285 t+1775 24416.24957+464.996s5+264.996 


s°+52s5’+926s5°+6700s°+208255°+25000s °+5611.06355" 


=-1 (6-74) 
-17555.746S- 17555.746 


or after factoring numerator and denominator yields: 





(s2+1.479s+1.519) (s+20.204) (s+4.316} (s+2) 


(s?+1.479s+1.519) (s+20.204) (s+4.316) (s+19.786) 


—————————————— (6-75) 
(st5.642) (s+1.4129) (s-0.840)’ 
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which is in agreement with the conclusions made in para- 
graph four, i.e., the number of the fixed roots is four 
which is also the order of the characteristic equation of 
the system before the singular cascade compensator was 
introduced and that the Root locus of the compensated 
Singular system with the parameter a used as the Root locus 
variable parameter has identical zeros and the same number 
of poles as the uncompensated system, i.e., the system 
before the singular cascade compensator was introduced. 

The Root locus of the compensated singular sys- 
tem was drawn by using a computer program in which equation 
(6-74) was implemented. This is shown in figure (6-31). 

In this figure the four fixed roots are shown with dark 


circles at the corresponding locations. 
a. Range of Dominancy of the Singular Root 


By the use of a DSL simulation program in 
which the system of figure (6-28b) was implemented, the 
range of dominancy of the singular root (S. = 0.60, 
oe  ™ 1.232586) was investigated. More than one hundred 
properly selected parameter pair values were used within 
ene ee Stable region. The conclusion from this 
analysis was that the singular root has a very wide range 
of dominancy. The detailed results which were interpreted 


on both, the s-plane and the a-f8 parameter plane, have as 


follows: 
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(1) S-plane. The Root locus of the com- 
pensated singular system is shown in figure (6-31). This 
Root locus has been drawn with the parameter oa used as the 
Root locus variable parameter. According to equation 
(6-72) and (6-75), the compensated singular system has 
eight characteristic roots, four of which are fixed at the 


points: 


S ~4.316 S ~0.7396+j0.9861 


1 


-20.204 S =02/7596-70. 93601 


ee 
1.e€., two fixed real roots and the fixed complex pair of 
the singular roots. All these fixed roots are shown in the 
corresponding locations of figure (6-31) with dark circle 
marks. The location of each of the other four character- 
istic roots is varying on the corresponding Root locus 
Segments as a function of the parameter a value. It must 
be noticed that a specific value of the parameter a defines 
a unique pair of parameter a and 8 values since for each 
Value of a corresponds one and only one value of 8 as can 
be seen from equation (6-73), which is the equation of the 
Operating singular line. | 

Interpreting the computer results which were 
obtained from a DSL simulation and a polynomial root finding 
program, the following remarks and conclusions can be 
Stated in correlation of the parameter a value and the 


Root locus of figure (6-31). 
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—= 7 


(a) O< a < 66.25. Within this range of 
the parameterq value, the compensated singular system is 
unstable because there is a real root on the right half of 
the s-plane. Based on the need for a specific name to be 
given to this root, Since it was found that after the 
Singular root this is the only root which affects the 
transient response of the compensated singular system and 
the fact that this root lies on the same Root locus segment 
with the "singular root", it was called the "non singular 


root’. 


(b) a = 66.25. For this value of the 
parameter a the non singular root lies on the origin of 
the s-plane and the compensated singular system is still 


unstable or conditionally stable. 


(c) 66.25 < a < 81.30. Within this range 
of the parameter a value, the non singular root lies on 
the real axis in the left half of the s-plane. The compen- 
sated singular system is stable and the non singular root 
is dominant. In figure (6-32) the, transient response of the 
compensated singular system is shown for a unit step input 
and with an arbitrarily selected value of the parameter a 
within this range, i.e., a = 75 for which the corresponding 
value of 8 was found from the equation of the operating 


Singular line, (equation (6-73)) to be 8 = 28.987. 
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(a) a 5-951. 350: For this value of the 
parameter a the non singular root is repeated at the 
point where the Root locus breaks away from the axis 
of reals, i.e., at this point the non singular root 
meets the other real root which was moving to the right 
on the real axis. This root is again the dominant root. 
In figure (6-33) the transient response of the compensa- 
ted singular system is shown for this value of the 
parameter a. As the value of the parameter a further 
increases, these two roots form a complex pair of roots. 
In this case by the name "non singular root" this complex 


pair will be meant. 


(e) 81.30 < a 95. Within this range 
of the parameter a value the non singular root is 
complex and although it has moved away from the real 
axis it is still closer to the origin and to the imaginary 
axis than the singular root. Within this range the con- 
tribution of the Singular root to the transient response 
of the compensated singular system was negligibly small, 
and the non singular root was dominant. In figure (6-34) 


the transient response of the system is shown for a = 90. 
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(£) 95 < a < 105. Within this range of 
the parameter qa value the non singular root has ap- 
proached closer to the singular root. The contribution 
of both of these roots to the transient response of the 
compensated singular system is significant and none of 


these can be said to be dominant. 


When the second order model cor- 


responding to the singular root (c. = 0.6,0,. = 1.2326) 


was considered, with transfer function: 


64) 
ns 1.5193 (6-76) 


2 2 2 
S#25.0, Stu. S05 739055295 


the corresponding transient response characteristics were 


calculated and they were found to be as follows: 


Time to peak of first overshoot: & = 3,19 sec 
settling time te 29.41 sec 
Percent overshoot : Mp = 9.48% 
Transient oscillatory frequence: w= 0.99 rad/sec 
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Within this range of parameter a value the transient 
response characteristics of the compensated singular system 


were found to be as follows: 


Time to peak of first overshoot: = 5 sec. 

Settling time : yas Eos oS ee 
Percent overshoot : 8 < Mp < 11% 
Transient oscillatory frequency: 0.6 < wy s 1 rad/sec 


Because in practical applications the 
transient response characteristics which are of the basic 
interest are the settling time and the percent overshoot 
it seem from the above results that an acceptable solution 
to the problem can be derived within this range of the 
parameter ao. In figure (6-35) the transient response of 


the system is shown for a = 100. 


(coerce howe se thin, Ghulss Gange oF 
the parameter a value the singular and the non singular 
root are very close to each other. The contribution of 
both of these roots to the transient response of the com- 
pensated singular system is again significant, and none 


of these can be called dominant. The transient response 
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was in general characterized by long settling time, which 
was approximately t, = 8.5 sec and high percent overshoot 
values, which were found to lie in the range from 11% up 

to 18.5%. In figure (6-36) the transient response of the 


compensated singular system is shown for a = 120. 


(h) a= 131.5. For this value of the 
parameter a the non singular root becomes identical with 
the singular root, i.e., the singular root is repeated. 
In figure (6-37) the transient response of the system is 


Shown for this value of the parameter a. 


(2) 2 ESie Se<-o-< 285." Within this tance 
of the parameter a value the non singular root, although 
it is further away from the origin and the imaginary axis 
than the singular root, its contribution to the transient 
response of the compensated singular system is still 
Significant and none of these roots can be called dominant. 
The transient response characteristics of the compensated 


Singular system were found to be as follows: 


Time to peak of first overshoot: 2.8 < bd < 4 sec 
Settling time : 7 < te < 8.8 sec 
Percent overshoot : 18.96 < Mp < 23.6% 
Transient oscillatory frequency: 0.98 < 4 < 1.1 rad/sec 
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In figure (6-38) the transient response of the compensated 


Singular system is shown for a = 200. 


(j) 285 < a < 385. Within this range of 
the parameter a value the non singular root has moved 
further from the imaginary axis and from the origin, and 
the singular root has become the dominant root but with 
low degree of dominancy. Within this range, the transient 
response characteristics of the compensated singular system 


were found to be as follows: 


Time to peak of first overshoot: ty 2G 

settling time : 6.6 < t, < 7 sec 
Percent overshoot : 14.22 < Mp < 18.96% 
Transient oscillatory frequency: 0.99 <¢ w, <1.1 rad/sec 


1.e., the percent overshoot is from 50% up to 100% over the 
corresponding to the singular root second order model 
value, and the settling time from 22% up to 30%. In figure 
(6-39) the transient response of the compensated singular 


System is shown for a = 330. 


Ghomeesopec oc <sloZ0.  Wathin this range 
of the parameter a value the non singular root is so far 
away from the origin in comparison with the singular root, 
that its contribution to the transient response of the 
compensated singular system is either very little or 


neglegible and the singular root is the dominant root. 
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Within this range of the parameter a value the transient 
response characteristics of the compensated singular sys- 


tem were found to be as follows: 


Time to peak of first overshoot: Vile) as ty Gon esec 
Settling time 9-6 < t, < 6.6 sec 
Percent overshoot 11.85 < Mp < 14.22% 
Transient oscillatory frequency: 0.98 < wy <0.99 rad/sec 


i1.e., the percent overshoot from 10% up to 50% over the 
second order model value and the settling time no more 

than 22% over that corresponding to the second order model 
value. In figure (6-40) through (6-44) the transient 
response of the compensated singular system is shown for 
several values of the parameter a within the range which 

was considered. It was noticed in this case that although 
the parameter values vary within a wide range, the transient 
response of the compensated singular system remains almost 


unaffected. 


(1) 1520 < a < 1922. Within this range 
of the parameter a value the non singular root appears again 
to affect appreciably the transient response of the com- 
pensated singular system due to the fact that although it 
is far away from the origin it is very close to the imagin- 
ary axis. The transient response is characterized with 
long settling time and high frequency of oscillations. In 


figure (6-45) through (6-47) the transient response of the 
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compensated singular system is shown for three selected 
values of the parameter a within the range which was 


considered. 


(m) oa = 1922. For this value of the 
parameter a the non singular root lies again on the 
imaginary axis and the system is unstable or conditionally 
Stable. The corresponding value of w for the non singular 


root at that point is w = 8.67 rad/sec. 


‘ 


(n} a> 1922. Within this range of the 
parameter a value the non singular root lies on the right 
half of the s-plane and the compensated singular system 
is unstable. In figure (6-48) the transient response of 
the system is shown for such a value of a,i.e., a=2000. 

The conclusion from the above analysis is 
that there is in general a wide range of the parameter a 
values for which the singular root is dominant, and which 
may be restricted according to the tolerance of the 
transient response characteristic values which are usually 


determined with problem specifications. 


(2) Parameter Plane.» In figure (6-49) the 
parameter plane diagram of the compensated singular system 
is shown for the singular values of Ce and Ws together 


with the singular line which is associated with them. As 
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can be seen in this figure all the above mentioned lines 
are tangent at the same point of the parameter plane. 

In figure (6-50) the same parameter plane 
diagram is shown for different scale values. On the same 
diagram the compensated singular system stable region and 
dominant singular root region have been interpreted in 


terms of singular line segments. 


b. Hardware implementation of the compensated 

Singular system 

According to Part C of this Section, the com- 
pensated singular system can be implemented as shown in 
figure (6-280) or (6-28b) or (6-28c), all of which are 
equivalent. 

Choice of the configuration to be used in 
practice will be in general determined from the specific 
application and the type of compensator available, i.e., 
integrated form or distinct components. 

When it is required to maintain the two para- 
meters a and 8 adjustable or when the ponerse ter of the 
system 1S going to be built from distinct components then 
the configuration of figure (6-28b) seems to be the most 
representative one for such a case. 

Considering also the advance of the technology 
today in integrated circuits, it must be realized that 
any required compensator can be built in integrated form. 


In such a case either figure (6-28a) or figure (6-28c) 
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can be implemented, after the two parameter values have 


been determined. 


c. Remarks 

When the singular line theory is considered as 
a design method for compensation of linear systems then 
the following remarks can be stated concerning the singular 
case which was considered in this section. 

(1) It can be applied to any type of plant. 

(2) The design problem has been simplified and 
better organized since in general the basic design problem 
is not to find the required compensator but the appropriate 
values of the two parameters a and 8B which have to be used 
with the singular compensator which was defined by equation 
(6-18). 

(3) The existence of the singular line pro- 
vides increased flexibility in the design procedure. 

It must also be noticed that when the two 
paramaters a and B are considered as time variant quantities 
of an adaptive system and its compensator respectively, then 
by adjusting their initial value and their rate of change 
in such a way that the operating point moves along the 
corresponding to the singular root dominant root region, 
System performance will be essentially unchanged. In 
this case other characteristic roots of the system will 
change but since they are not dominant, system performance 


can not be affected. When the above mentioned rate 
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adjustments are made automatically the compensated singular 
system becomes self-adaptive. This type of application 

of the singular line theory in adaptive systems needs to 

be further investigated. 

In the previous example the results 
which were obtained were related with a mathematical 
accuracy which in general can not be achieved in practice. 
In this case a reasonable question will be: How is the 
performance of the compensated singular system affected 
when the values of the two parameters a and 8B are set 
approximated with respect to their theoretical values? 
or what is the sensitivity of the singular dominant root 
with respect to each or both of the two parameter values? 


This problem is next analyzed. 
G. SENSITIVITY ANALYSIS 


if Introquction 


The characteristics of a physical component are 
Subject to change for various reasons. Perfect control 
systems with stable components which can be built or can 
maintain their design values under any conditions are 
eee to obtain in practice. Take for example such 
a passive component as a resistor. It is almost impossible 
to obtain a resistor with the exact design resistance 


value. For this reason a tolerance usually ranging from 
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7% up to 30% will accompany its ‘specifications. Even if 
this "ideal'' component is obtainable, factors like tem- 
perature variations, radiation etc., may eventually cause 
the resistance value to vary and produce undesirable 
results on the control system performance. These general 
rules are also applied when an active or solid state de- 
Vice is considered. 

A measure of the effects of parameter variation 
upon the performance of a control system has in general 
been demoted as sensitivity. Basically there are two 
typesof parameter variations studied by various inves- 
tigators. The first one which is the most popular and 
the most common in practice is the incremental variation 
case, where the change in parameter values is assumed to 
be very small, i.e., infinite-simal. The corresponding 
sensitivity study in this case is called <<Microscopic 
Sensitivity Analysis >>. The second case is that of 
relative large variations in parameter values l.e., 
finite changes. In this case the corresponding sensitivity 
study is called << Microscopic Sensitivity Analysis >>. 

The definition also of the sensitivity coefficient 
varies among investigators. Generally the sensitivity 
coefficient S can be defined as the ratio of the amount by 
which a control system's performance characteristic deviates 
from its nominal (or original) value when one or more sys- 


tem parameters vary to its nominal value. From this 
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definition the sensitivity coefficient is mathematically 


given by: 


gs 2 Ac (6-77) 


where C denotes the nominal value of the system charac- 
teristic which is studied and Ac the deviation from this 
value when one or more system parameters which are con- 
Sidered vary. 

There have been many sensitivity functions 
defined for the effects of parameter variations on the 
transfer functions of linear control systems, some of 
which are the following: 

a. The Kokotovic Method 

The root sensitivity functions developed by 
Kokotovic and Siljak provide most general solution to the 
sensitivity problem for small parameter variations in 
linear control systems and the most applicable to parameter 
plane methods. 

b. The F.H. Hollister Method 

The macroscopic root sensitivity functions 
defined by F.H. Hollister are applicable for large para- 
meter variations. 

c. Bode Method 

Bode defines a logarithmic or normalized 
sensitivity function of a variable, Xs Gt: Se with 


respect to a parameter 25% asi. 
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iL 
In general an analytic investigation of the 
sensitivity problem of a linear control system is neither 
Simple nor easy. By computing sensitivity coefficients 


limits of individual parameter variations may be deduced 


and tolerance levels can be formulated. 


2. Sensitivity Analysis of a Singular System 


The sensitivity analysis of a singular system 
can be performed by any method applicable to a linear 
feedback control system. By such a method the sensitivity 
of any system performance characteristic due to variations 
of one or more parameter values can be analyzed. Especially 
when compensation of a control system is designed by the 
use of the "singular line theory", then the transient 
response characteristics and the two variable parameters a 
and 8 are of the main interest flor the singular system 
which results, at least for the first steps of the design 
process. Therefore it will be very important to know how 
transient response characteristics will be affected and 
how tolerance levels can be established concerning varia- 
tions of the two parameter a and 8 values. In terms of 
sensitivity this could be stated as: How. sensitive are 


the transient response characteristics of a singular system 
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to variations of the two parameter values and how can the 
corresponding tolerance levels be determined? 

The goal of the present analysis is to give an 
answer to these questions by a relatively simple and 
accurate way, in which computer use will be optional. 

Generally when designing by the use of 
Singular line theory, specifications on transient response 
characteristics determine the desired operating singular 
point (see paragraph VI, F, 6). From this singular point 
or singular root the equation of the corresponding singular 
line will be evaluated. Some other specifications or 
physical limitations of the system may also restrict the 
permissible operating region on the singular line to a 
finite line segment. Then when the operating point moves 
or stays fixed on this singular line segment, performance 
requirements on transient response characteristics are 
satisfied. The design process up to this point has been 
described in detail in the previous part of this section. 

It will be assumed now, that the operating point 
moves away from the singular line due to variations of 
either one or both of a and 8 parameter values from their 
design values. This will be in general the case for an 
actual system. In such a case the acceptable tolerance 
on the transient response characteristics will determine 
how far away from the singular line the operating point 


1S permitted to move. 
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Sincemtiememloeat least one singular Pine for 
each value of C for a singular system it will be expected 
that for a singular system there is an infinite number of 
Singular lines on the parameter plane, each of which 
corresponds to a different value of %. Therefore accord- 
ing to the tolerances, which are given for the transient 
response characteristics of the compensated singular 
System, a permissible region can be established on the 
parameter plane for the operating point, which will be 
bounded between two such singular lines. This concept 
not only was applied in the following example but was 
also further investigated from a ieore rica point of view 
and it was found that it is true and valid, which will be 
shown later in the example. 

Determination of these two singular lines on the 
parameter plane will not only give the answer to the pre- 
vious question, but it will give also a lot of other 
information concerning sensitivity and performance of a 
Singular system for parameter variations. It will also 
give a precise way for tolerance values establishment for 
a singular control system. 

Although different specifications or different 
Systems may give rise to different boundaries for the 
operating aeeint on the parameter plane or different cri- 
teria in order to determine these boundaries, the procedure 


which must be followed will be in general the same. This 


will be illustrated by the following example. 
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Se Example | (6-7) 
Consider again the plant of example (6-6) with 
transfer function: 


20 


G(s) = ——--—-—_——_ 
S(st+1)(s+5) 


which 1S compensated by an initial cascade compensator 
with transfer function: 


s+Z 


G(s) = 
< s+2Z0 


The Root locus of the resulting system which is 
drawn with the gain of the plant used, as the Root locus 
variable parameter, is shown in figure (6-51). 

It is already known (see paragraph VI, F, 6) 
that any point of the above Root locus, which does not 
lie on the real axis is a singular point for the compen- 
sated singular system. It will be assumed that the 
specifications on the transient response characteristics 
for the compensated singular system require the operating 
point to be the Root locus point which is located at 
(. = 0.60, ean 1.23) and that according to the given 
tolerance on ene transient response characteristics which 
are interpreted in terms of the second order model, this 
Singular point is permitted to move within the finite Root 
locus segment which is bounded between the Root locus 


points: 
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! ! a 
4! 0.65, wi = 1.03) 


and 


CB. = 555, Whe ee lacy 

These points can be determined by the intersection 
of the Root locus with the corresponding constant ¢ lines 
(see figure (6-51)) or can be evaluated by using the 
"Singular point program" in the computer. In the second 
case, although more precise values can be obtained for 


u) t 


me and Was. it will be seen later in this example that 


Such an accuracy is not actually needed. 

Since the above ae eed Root locus segment has 
been established by interpreting problem specifications 
in terms of second order model, will be actually the 
permissible region for the operating singular root, if 
and only if, this singular root has a high degree of 
dominancy anywhere within this region. 

By, either studying the Root locus of figure 
(6-51) or taking into account the results obtained from 
example (6-6), a basic remark concerning the dominancy of 
the singular root in the case which is considered is that 
as 6. decreases, the degree of dominancy of the singular 
root also decreases because it moves away from the origin 
and from the imaginary axis, and the contribution of the 


non-singular root (see example (6-6)) becomes more 
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Significant. Based on this remark and the need for the 
given specifications on the transient response charac- 
teristics to be met, it was decided the permissible 

region for the operating singular root locus to be restric- 


ted between the points: 


coe OOO, me 1.23) 


and 


(ee Wa ae We = 1.03) 


which are shown in figure (6-51). 

At this point it must be noticed that such a 
decision is based absolutely on personal judgments the 
correctness of which can be checked by simulation of the 
meipensated Singular system at the limiting acceptable 
values for the two parameters a and 8, which were obtained 
from this analysis, as it will be shown later in this 
example. 

By introducing a singular cascade compensator, 
the characteristic equation of the compensated singular 


system becomes: 


ee orice (GOO ans Maze 2Se2en) s+ 


(25000+177a)s°+(100004+3500+208)s*+ (6-78) 
(200a+808)s + 808 = 0 
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When the singular point (oc. = 0.60,w ). = 1.23) 
was considered, the equation of the corresponding singular 


line was found from equation (2-5) and (2-16) to be: 
B = 3.3040 - 218.355 (6-79) 


Comparison of, equations (6-79) and (6-73) reveals 
that the use of the computer for evaluation of the polar 
coordinates of a singular point and from these the equation 
of the corresponding singular line to be calculated, does 
not give any significant improvement a the accuracy of 
the equation of the singular line, for analysis and design 
purposes. 

When the singular root Ca = 9.65, We = 1.03) 
was considered the equation of the corresponding singular 


line was found from equation (2-5) and (2-16) to be: 
B = 2.4444 - 138.097 (6-80) 


The above two singular lines were drawn on the 
parameter plane diagram which is shown in figure (6-52). 
After taking into account the results of example (6-6), the 
Mracteristic equation of the compensated singular system, 
and equations (6-79) and (6-80) of the boundary singular 
lines, the stable region ABCA BCA, and the region 
DE F G in which the varying singular root has a high 
degree of dominancy were approximately established on the 


parameter plane diagram of figure (6-52). 
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When the operating point lies at the point A 
where the two singular lines intersect, then the compen- 
sated singular system has both the corresponding Singular 
roots as characteristic roots. When the operating point 
lies within the region D E F G where the varying singular 
root is dominant with a high degree of dominancy, then 
the specifications on the transient response characteris- 


tics for the compensated singular system are met, i.e: 


Peak overshoot Mp: From 6.8% up to 12.6% 
settling time ty: From 4.6 sec up to 5.9 sec 


Trans. oscil. freq. we? From 0.78 rad/sec up to 1.32rad/sec 


This was actually verified in the computer by a 
DSL simulation program of the compensated singular system, 
in which several pairs of a and 8 parameter values which 
define points in the region D E F G of the parameter plane 
diagram of figure (6-52) were considered. In figure (6-53) 
through (6-56) the time response of the system for some 


of these points is shown. 
a. Sensitivity Study 


Each operating point within the region DEF G 
of figure (6-52) lies on a singular line (it will be shown 
later in this example), which corresponds to a singular 


root (5520, ). These polar coordinates of this singular 


S 


root can be approximated for any such point from the 
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corresponding values of the boundary singular lines. For 
example the point M, and any other point which lies half 
way between the boundary singular lines will correspond 
to a singular root (o. =<On02o% Wis - 1.13). When such a 
Moint, i1.e., the point M, (ao = 900, B = 2400) is 
considered the corresponding transient response charac- 
teristics can be evaluated from the second order model 
which corresponds to the singular root (5. i 0256 Wis - 
13). 

Considering now that due to variations of 
elther one or both of the parameter a and 8 values the 
Operating point moves to another location within the 
region D E F G of the parameter plane of figure (6-52), 
then the corresponding new G and Wg Values can be 
approximated and from these the corresponding transient 
response characteristics. Therefore the sensitivity 
coefficient which relates the transient response 
characteristics and the variations of the two parameter 
a and £6 values can be evaluated and from these sensitivity 
Studies can be made as it is shown below: 

(1) Sensitivity of the Peak Overshoot due to 

Variations of the Parameter a. Let the 
operating point move from M, Co = 900, 6 = 2400) to M, 
(ao = 990, B= 2400) due to a +10% change of the parameter a 
value. The corresponding singular points were approximated 


as (o. = 0.625, Ones 1.13) and (52 =O Ole Ws so On) 
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respectively. The corresponding percent peak overshoots 
as they can be found from the second order model are 

Mp = 8% and Mp = 7.3%. Therefore according to the defi- 
nition of the sensitivity coefficient which was given in 


paragraph one: 


MP _ Mp'- 
_Mp'- Mp _ 7.3-8 _ _ 


(Zy een STELV ity On the Peak Overshoot due 
to Variations of the Parameter 8. Let 

now the operating point more from M, ( a = 900, B = 2400) 
to Mz (a = 900, 8 = 2640) due to a +10% change of the para- 
meter 8 value. The corresponding singular points were 
approximated as (o. = 0.625, Ws 1.13) and (oo = 0.61, 
Was = 1.19) respectively. The corresponding percent 
overshoot was evaluated from the second order model to be 
Mp = 8% and Mp" = 8.9%. Therefore the sensitivity co- 
efficient which relates the peak overshoot and the +10% 


variation of the parameter 8 value from its value at the 


mittial point M, will be: 


MP Mp - 
S ay pip Seo. = 0.1125 
B Mp 8 
Comparison of the above evaluated two sensi- 
tivity coefficients reveals that when the operating point 
M, is considered, then the peak overshoot of the compensated 


Singular system is more sensitive to variations of the 
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parameter § than it is to variations of the parameter a 
(almost 1.28 times more sensitive), and that the effect 

of the variations of these two parameters on the Peak over- 
shoot of the compensated singular system is opposite. 

It must be noticed that although the values 
obtained for the sensitivity coefficients are approximate, 
the conclusions which can be derived from a relative com- 
parison must be considered as very accurate for analysis 
and design purposes. 

In the same way any sensitivity problem which 
relates the transient response characteristics of the com- 
pensated singular system and variations of either one or 


both of the parameter a and 8 values can be studied. 


b. Selection of the Operating Point 


In general the operating point M of a compen- 
sated singular system must be selected where the Singular 
root has the highest degree of dominancy (which results in 
the best possible approximation of the transient response 
characteristics of the compensated singular system by the 
corresponding second order model) and on a point halfway 
between the boundary singular lines (which will give the 
maximum tolerance for variations of the two parameters a 
and 8, when it is considered that both can vary simul- 
taneously). Under these considerations the best location 


for the operating point in the case which was studied in 


134 





this example seems to be at the point M(a = 1310 , 8 = 
3500), which is shown in figure (6-52), It must be noticed 
that some other specifications of system's performance, or 
other limitations may require another location for the 


operating point M to be determined within the region 


mE FF G. 


c. Tolerance Levels Establishment on Parameter a 


and 8 Variations 


After the operating point M has been chosen, 
then tolerance levels on the two parameters a and B 
variations must be established. These tolerance levels 
must be determined in such a way that the operating point 
will be restricted to lie within a circle which has center 
the point M itself and is tangent to the two boundary 
Singular lines and to any other boundary line in the 
neighborhood of this point. For the case of the example 
which was studied the required tolerance levels for both 
the a and 8 parameters are +7.5% from the corresponding 
nominal values of a = 1310 and 8 = 3500, as it can be 
evaluated from figure (6-52), i.e., for no more than 7.5% 
Variations towards both directions of each parameter a 
and 8 value, the operating point will lie within the circle 
Shown in figure (6-52) whose center is at the point M, 
and the performance of the compensated singular system will 


be such that the specifications will be met. 
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4, Comments 


In establishing tolerance levels in practice some 
trade offs can be made by. considering each parameter 
Separately since the nature of each is different, i.e., the 
parameter 8 is always associated with a gain value while 
the parameter a is associated either with another gain 
value or with pole locations (see figure (6-28a,b,c)). If 
for example, in the previous case, which was studied above, 
the poles of the singular compensator can be set very 
accurately, due to the nature of the singular compensator 
which 1s going to be used, then its gain 8 will not be 
required to be set with a higher than +15% accuracy, in 
order for the operating point to be restricted within the 
desired circle which was mentioned above, as can be seen 
from figure (6-52). 

Another general remark which was derived from 
figure (6-52) is that as far away from the intersection of 
the two boundary singular lines the operating point M 
1s located, the less sensitive the compensated singular 
System becomes to the parameter a and 8 variations. 

The above comments together with the results 
obtained up to this point from the sensitivity analysis 
reveal that design of compensation by using singular line 
theory leads to very accurate setting of tolerance levels 
for the two parameter a and 8 values and very easy and 


accurate way of determining the optimal operating point for 
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the compensated singular system either on the s-plane or 
on the a-8 plane. 

Although all the required information in order to 
make a conclusion about the sensitivity of a compensated 
Singular system transient response characteristics, due to 
Variations of the parameter a and 8 values, have not been 
derived, based on the results of the analysis which was 
performed up to this One it is believed that in general 
compensation of a system by using singular line theory 
may lead to control systems which are not sensitive or at 
least very sensitive to parameter values variations. 

In the previous example it was shown that when 
the operating point moves within the region DE F G of 
figure (6-52) then the transient response characteristics 
of the compensated singular system can be well approximated 
by the corresponding second order model. Actually there is 
no question about it, when the operating point lies on a 
Singular line, like the two boundary singular lines D E 
and G F of the previous mentioned operating region. In the 
case where the operating point lies between these boundaries, 
although the results from example (6-7) reveal that the 
method is still valid, some reasonable questions which may 
rise, in this case from a theoretical point of view are: 

Is the compensated system still singular? Do singular 
conditions exist? Does the operating noi Stsull lie on 


a singular line? or have all the characteristics of the 
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Singular conditions been lost? These questions were 
answered as follows: 

When a linear control system with negative unity 
feedback, whose open loop transfer function is G(s) = 
ES is compensated by a singular cascade compensator, 
the characteristic equation of the compensated singular 


System will be: 
1 +a F_(s) +8 KF a) = 0 
& & 


It has been shown (see paragraph III, C, 3) that such a 
System has in general an infinite number of singular points, 


1.e., points of the s-plane which satisfy equation: 


/F,(s) = nil 


where n is an real integer or zero. When such a singular 
point has been specified, the equation of the corresponding 
Singular line can be determined either by direct substitu- 
tion of the value of the singular point s in the character- 
istic equation of the compensated singular system or by 
using equations (2-5) and (2-16). From the equation of 
this singular line then, an infinite number of a and 8 

pair values corresponding to the singular point under 
consideration, can be determined. In the opposite case 


Where a pair of ao and 8 values has been specified the 
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equation of the corresponding singular line and correspon- 
ding to this equation the singular point can also be 
determined. 

This will be illustrated by considering the follow- 
ing two cases from example (6-7) where the operating point 
moves within the region D E F G. In the first case the 
operating point arbitrarily was assumed to lie at the point 
M,(a=900, 8=2400) of figure (C_ S72 euetiren. the Characteristic 


equation (6-78) of the compensated singular system becomes: 


6 4 5 Z 


s°+525/+926s +76008°+46025s +184300s° +373000s 


+372000s+192000 = 0 


By the use of a computer program with the ZPOLR 


subroutine the roots of the above equation were found to 


be: 

Sy = - 0.7044558+j0.8843676 S5 = - 0.7044558-j30.8843676 
a > 0.8651241+j55.8923703 eo = 0.8651241-45.8923703 
Se = - 22.1206845 S¢ = - 20.1762447 

So = - 4,4148436 S = - 2.1490673 


Each of the above complex roots was checked and 
it was found that each lies on the Root locus of figure 
(6-51). Therefore each of these complex roots is also a 


Singular root, (see statement after equation (6-61)). A 
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more general theoretical proof of this statement which 

is always true, can be derived from paragraph VI, F, 4. 
Because the operating region of the parameter 

plane, where the operating point My lies has been designed 

as dominant root region only the first singular root, i.e., 

the singular root Sy9> which also was dominant, was further 


considered. This root is: 


S, > = - 0.7044558 +} 0.8843676 


or 


(ae Mic Sreao eos: Talo) 


When this singular root was considered, from 
equation (2-5) and (2-16), the equation of the corresponding 


Singular line was found to be: 


Bee c O42 be LOZ. 1555 (6-81) 


Since the coordinates of M,> i1.e., a = 900 and 
B = 2400 satisfy the above equation it implies that this 
point lies on a singular line of a dominant singular root 
whose coordinates are (S. = 0.623, ee ee Sie) 

Two points are interesting here. First, the above 
coordinates of the dominant singular root were approximated 
from the corresponding values of the two boundary singular 
lines of figure (6-52), during the sensitivity study of 


example (6-7), (see paragraph G3 a (1) of this section), 
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and second, when the other singular root was considered, 
then by the same method, it was found that the operating 


point M, 


this root, i.e., finally the operating point My lies on 


the intersection of two singular lines. Of course under 


lies also on the singular line corresponding to 


dominant root considerations this last singular line is 
of no interest. 

It was next assumed that the operating point was 
moved to a new location M; Co = 1430,8 = 4280) of the 
dominant root region D E F G of figure (6-52). Then the 
characteristic equation of the compensated singular 
system was evaluated to be: 

594525 /+9265°4+8130s°+60825s 42781108 +596100S~ 
+628400 s + 342400 = 0 


] : 
By the same method as in the previous case the roots 


of this equation were found to be: 


= - 0.7257602+50.9470602 Ss, = - 0.7257602-j0.9470602 
Sz = - 0.3981063+j7.5192493 By, 6) = Wo cE UNOO S31! SIGs 
eee - 23.1172585 s. = - 4.3550711 
s, = - 20.1934085 Sg = - 2.0865288 


Each of the above complex roots was also checked 


and it was found that it lies on the Root locus of 
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figure (6-51). Therefore they are also singular roots. 


In this case the dominant singular root is: 


S$, 9=- 0.7257602 +j30.9470602 
? 


or 


Ge 2 O08, mag = ones) 


\ 


When this singular root was considered, from 
equation (2-5) and (2-16), the equation of the corresponding 


Singular line was found to be: 


B = 3.13070 - 196.1499 (6-82) 


Since the coordinates of the operating point Me 
which was considered, i.e., a = 1430 and f= 4280, satisfy 
the above equation, it implies that the operating point 
Me lies on the singular line which is defined by equation 
(6-82). It must also be noticed that equivalent remarks 
like these which were stated in the previous case about the 
Other complex root 334 also were found to apply and in 
this case. 

By the same way the singular lines and singular 
points corresponding to any point within the region DE FG 
Or to any other point of the parameter plane, can be 
evaluated. 


The conclusion from these results is that indepen- 
dently of where the operating point lies on the a-8 plane, 


Singular conditions always exist. 
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The results from the example which was studied are 


summarized in Table I of Appendix Ga 


H. STEADY STATE ACCURACY OF A SINGULAR SYSTEM 


ine introduction 


A measure of control system accuracy is usually 
the steady state error. Related concepts and definitions 
can be found in reference [2]. 

When steady state error considerations are made 
for a singular compensated system shown in figure (6-5), 
the question which arises is: What is the effect of the 
Singular cascade compensator on the error constant LS 
(which is also called the DC gain or zero-frequency gain), 
of the system which is compensated by a singular cascade 
compensator? 

The goal of the present study is for an answer to 
be given to this question and for the required design 


procedure to be established. 
2. Steady State Error Analysis of a Singular System 


It is known that if a finite error exists in steady 
state operation of a control system its magnitude is deter- 


mined by the reciprocal of the error constant Ky: 
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When a unity feedback control system is considered 


with open loop transfer function: 


(s+z )(s+z ) vee (S+Z, ) 
G(s) = k F,(s) = k a 
s (stp J(stp,) .-. (stp) - 


then the error constant K 1S: 


K. = gim s G(s) (6-83) 
S20 
or 
ee ee (6-84) 
x ‘ 
P,P, + ++Py 


When the above considered system iS compensated 
by a singular cascade compensator whose transfer function 
G..(S) iS given by equation (6-18), then the error constant 


of the compensated singular system is: 


22 ae N 7 
ae = Lim s G..(s) G(s) (6-85) 
S>O 
or 
kB zz Z 
= Tae? eee ean : 
eee : : - KY (6-86) 
ee Pa 
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i1.e., the error constant of the compensated singular 

system differs from the error constant of the system 

before the singular cascade compensator has been introduced 
by a factor equal to the ratio 8/a of the two parameter 
values. Since an operating singular line does not in 
general pass through the origin of the a-8 plane, the 

ratio 8/a will not in general remain constant when the 
operating point is moved along a singular line. 

When an operating region on the a-§8 plane is 
considered, which actually will be the case, i.e., the 
parameter values will be fixed, then the error constant 
of the compensated singular system will have also a fixed 
value. 

Since the ratio §8/a corresponds to the slope of a 
straight line through the origin of the a-8 plane, it 
implies that the locus of all the points of the a-§8 plane 
with a fixed value of the ratio 8/a is a straight line 
through the origin, and that the locus of all the points 
of the a-8 plane with a specified range of values of the 
ratio 8/a is a region of the a-f8 plane which is bounded 
between two such straight lines through the origin. 

Therefore interpretation of any steady state 
accuracy specifications in terms of the 8/a ratio value 
leads to "steady state error boundary lines" establishment 


on the a-g plane as is shown in the following example. 
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3. Example (6-8) 

In this example it was assumed that steady state 
accuracy specifications were established for the compensated 
Singular system of the previous example (6-7), and that 
interpretation of these in terms of the error constant 
K set the following restriction on the 8/a ratio value, 


SX 


1.e.: 


—£ >] (6-87) 


The corresponding "steady state error boundary 
line’ for this case is the straight line through the origin 
with slope equals to one. This was drawn on the a-8 plane 
shown in figure (6-52). Then all the points of the a-8 
plane above this line satisfy steady state accuracy 
Specifications of the problem. Therefore any point of the 
dominant singular root region D E F G shown in figure (6-52) 
will also give an acceptable steady state error constant 
for the compensated singular system. 

In the same way, when another system or different 
Specifications are given concerning the steady state 
accuracy establishment of the corresponding "steady state 
error boundary lines" or line, on the parameter plane the 
required information for analysis and design purpose can 


be derived. 
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Vad; THIRD SINGULAR CASE ANALYSIS 


In this section the "third singular case" which was 
defined in Section IV, was considered for analysis. This 
Singular case corresponds to the compensation structure of 
figure (4-5) where: 

G(s) = st) 


Gc(s) = ps F, (s) (7-1) 
H(s) = as 
and the plant G(s) = kF g(s) has at least two poles at 


the origin, i.e., the cascade compensator has no zero's 
at the origin. 
The characteristic equation of the compensated 


Singular system is: 


ltaks F,(s)+8k (sF,(s)4) = 0 
A. SYSTEM OUTPUT AND ERROR ANALYSIS 


The transfer function of the compensated singular 


system is: 


-<2R 2 
T(s)= C(s) _ Bks Fo (s) 


R(s) 


72) 
l+aks F (s)+8ks“ BoA {s) 
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and the output and. error are respectively; 


C(s) = R(s) TCs) Ci). 
and 
E(s) = 


R(s) {1 - TC(s)} (7-4) 


1. Plant Type Two 


When the plant is type two, i.e.: 


(s+z,) (stz,)...(stZ,) 


G(s) = kFg(s) = k—~——___*+_____*__ (7-5) 
s'(st+p,) (s+p,)..-(s*p,) 
equations (7-3) and (7-4) become respectively: 
m 
ek [ ]cs+3)- 
C(s) = Se 7 


Nl 


Mm qh ™m 
s[ | (stm teks [ Ics) F Vcsemdeak [ Ies+25)° 
k j=1 k=1 j=1 


(7-6) 
and 
, Tl ‘ m Nn 
s[_[(sty) teks | lest) | lesen 
E(s) = re ag 
S | |cs+p,) +aKs | |(s#zg) | | (s+pK) +8k Ga) 
k=1 j=1 k=] j=l 
Cai?) 


When the final value theorem was applied to the 
above equations the steady state output and error were 
evaluated assuming that all the poles of the transfer func- 


tion T(s) lie on the left half of the s-plane. 
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a. Step Input 


When a step input R(t) = Au(t) or R(s) = A/s 


was considered it was found that: 


C(t) = Ceeiee) = A 
tro 

and 
E(t) = E,,(t) = Q 
ta 


b. Ramp Input 


When a ramp input R(t) = Btu(t) or R(s) = B/s¢ 
was considered it was found that: 


Gee) = Ghee Cele = 


t+« 
and 
| 
Bas, . 'p 
E(t) = E..(t) = seu 
tro SS m 
a! | z. 
k= 
c. Parabolic Input 
When a parabolic input R(t) = Ct7u(t) or 
R(s) = NG was considered it was found that: 
iG ie Ce CE ne 
t7o 
and 
E(t) = E.,(t) = « 
ee 
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2. Plant Type Three 


In the same way, the steady state output and error 


value for a type three plant were found to be as follows: 


a. Step Input 


Oe) ae) aK 
oe 3 
and 
Gee E Ct) = 0 
ter Ses) 
b. Ramp Input 
C(t) = C__(t) #e 
t+ = $5 
and 
n 
Bal |p 
Cee i ae ee ae 
Ss m 
t+ = 
al |e, 
ina 
c. Parabolic Input 
CEt) = C.,(t) = 106 
t7a 
and 


pe = E(t) Tlie « 
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3. Remarks 


The general form of the characteristic equation 
of the compensated singular system in the case which was 


studied in this section is: 


2 
ltaks FCs) + Bkis 5 (s)} = 0 (7-8) 


where the transfer function of the plant is: 


i 


G(s) 
By defining: 
Gs) a eshs (Ss) 


the characteristic equation (7-8) becomes: 


KE, (s) (7-9) 


Meesiaki(S eb lC UE - (ai =a (7-10) 
which is the characteristic equation of a singular system. 

Comparison of equations (7-10) and (6-16) reveals 
that the design process which has to be followed when 
compensation is performed by the compensation scheme, 
which corresponds to the "third singular case" will be in 
general the same with that corresponding to the "second 
Singular case''. The only difference will be that instead 
of using the Root locus of the uncompensated system as the 
basic design tool, the Root locus of the uncompensated 
System with one pole less than the actual number of poles 
at the origin has to be used. This will be illustrated in 


the following example. 


1S 





B EXAMPLE (7-1) 


In this example the plant shown in figure (7-1) was 
considered with open loop transfer function: 


1 
s (st2) 


1 


G(s) = KF, (s) 100 (C7151) 
This plant is unstable and it was assumed that prob- 
lem specifications dictate for a stable system with a 
dominant root with polar coordinates: 
c = 0.40 


and 


2<w <3 


The ponnenenred Singular system is shown in figure 
(7-2). As can be seen the whole design problem consists 
of finding the appropriate pair of ao and g parameter 
values by which the given specifications can be met. This 


was performed in three steps which follow: 


ie ocep, one 


f 


thie fins: Step Was to determine the required 
operating singular root for the compensated singular sys- 
elm. 

The Root locus of the uncompensated system with a 


pole less than the actual number of poles at the origin 


I Bye 





is shown in figure (7-3), i.e., this is the Root locus of 


equation: 


ae ee ae) 


5((SNN 


By the same reasoning as in the case of the "second 
singular case" it can be shown that any point of the above 
Root locus, which does not lie on the real axis is a sin- 
gular point of the compensated singular system. (See Section 
VI, F, 4, b). Considering next the specifications the 
desired singular point was located on the Root locus at the 


point (c = 0.40, es tee 


2. Step Two 


The second step was to find the equation of the 
Singular line, corresponding to the selected singular 
point. 

From equation (7-9) the characteristic equation 


of the compensated singular system was found to be: 


Be cA OMene 20s aloe =o Gas) 


When the singular root UE A 0.40, w = 2.5) of the 


ns 
above characteristic equation was considered, the equation 
of the corresponding singular line was found from equation 


(2-5) and (2-16) to be; 


B = 6.250-0.390625 Cau 
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-_ 


3. Step Three 

The third step which was the last one was to find 
the required operating point of the singular line which 
was defined by equation (7-14). 

Because it was also desired that the design pro- 
cess for this problem be completed without using a computer 
this step was performed as follows. 

Introducing equation (7-14) into (7-13) the 
characteristic equation of the compensated singular system 


becomes: 


Bae AA eOUa is - H200cS+6250.59.0625 = 0 7s) 


Although according to the results which were stated 
in Section VI, F, 4, the Root locus of the compensated 
Singular system, i.e., the Root locus of equation (7-15), 
with the parameter a used as the Root locus variable para- 
meter has the shape of the Root locus of equation (7-12) 
shown in figure (7-3), and therefore the drawing of the 
Root locus of equation (7-15) is not actually required, it 
was done in this example for a more detailed illustration. 

Equation (7-15) can be written as: 


100s“+200s+625 Sel 
jie eS eee 


5445 +45°-39.0625 
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or after factoring numerator and denominator; 


5 ieCaicn cob) jes Cl- 12.290) } ae 
{s+(1+j2.291) } {s+(1-j2.291)} {s+3.693} {s-1.693} 


The Root locus of the above equation is shown in 
figure (7-4). This root locus which was drawn by hand, 
consists of a fixed point there where the singular root 
(7 ,=0.40, Wig nt: 9) or s#-1+j2.291 is located and the 
regular Root locus segments. 

Comparison of figures (7-3) and (7-4) reveals that 
the two Root loci have actually the same shape. 

Because the singular root and the complex non- 
Singular root (See Section VI, F,4,a,(1),(c)), have always 
equal real parts, in this case, in order for the singular 
root to be dominant its residue must be much greater than 
that corresponding to the non-singular root, i.e., the 
Singular root must lie much closer to the origin of the 
S-plane than the non-singular root. Therefore the non- 
Singular root must lie far away above the singular root 
which is shown on the Root locus either of figure (7-3) or 
ie -4). 

By arbitrarily selecting the distance of the non- 
Singular root from the origin to be relative very large, 


in comparison with that which corresponds to the singular 
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root, dominance of the singular root in general shall be 
expected, e.g., let the non-singular root be: 
s = -1 + j50 
Then the left part of equation (7-15) must be equal to the 
product of the singular root and non-singular root. This 
yields a = 25.07. From equation (7-14) the corresponding 
value of the parameter 8, was found to be Bg = 156.31. 
Therefore the required parameter pair values or the 
desired operating point of the a-8 plane is (a=25.07, 


B=156.31). 


4. Remarks 


The compensated singular system shown in figure 
(7-2) was simulated in the computer by using a DSL program 
in which the two parameters a and 8 were set at the above 
indicated values. 

From this program it was found that the time 
response characteristics of the compensated singular system 
were exactly the same as these corresponding to the singular 
root second order model. It is shown in figure (7-5). 

It must also be noticed that the desired perform- 
ance of the compensated singular system can be achieved by 
an infinite number of a and 8 parameter pair values, which 
Berrespond to a different selection in the location of the 


non-singular root on the s-plane. 
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Criterion for the selection of the appropriate 
location for the non-singular root can be based either on 
sensitivity or and steady state accuracy considerations as 


follows: 


a. Sensitivity Considerations 


As the location of the non-singular root moves 
away from the axis of reals on the Root locus of figure (7-4) 
the value of the Root locus variable parameter a increases. 

Considering next the parameter a-f8 plane and 
the singular line (equation (7-14)), which corresponds to 
the operating singular root (z. Ses Olas 2.5), as the 
value of the parameter a increases the operating point 
(a,B) moves away from the origin on the singular line, but 
according to the results obtained from the sensitivity 
study in the previous section it yields a less sensitive 
compensated singular system with respect to variation on 


the parameter a and 8 values. These remarks can be used 


accordingly in each case. 


b. Steaily State Accuracy Considerations 


When the error constant of the compensated 


Singular’ system was considered it was found to be: 





Kee se 
Re 20 


When problem specifications concerning the 


error constant K , are also established, then an additional 


ILS) 








criterion for the selection of the location of the non- 
Singular root will be the desired value of the parameter «a 
for which the ratio 8/a is such that the optimal error 
constant coefficient can be achieved. 

In general the ratio 8/a can be effected either 
by selecting a different location for the non-singular root 
or by selecting a different operating singular point, i.e., 
operating on a different singular line or by introducing 
additional gain adjustments at the forward path of the 
minor loop of figure (7-2). 

[It must be noticed also that any other specifica- 
tions can be taken under consideration in each design step 


of the previous illustrated example. 
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Va. EXTENSIONS OF THE SECOND SINGULAR 


CASE COMPENSATION STRUCTURE 
A. INTRODUCTION 


Among the three singular cases, which were previously 
analyzed, the most interesting for applications in design 
of compensation of a linear feedback system by the 
dominant mode method, appears to be the "'second singular 
case" (See Section VI). 

The basic advantages of this case are that the plant 
type number is preserved, it can be used for any type of 
plant and for any plant structure. 

On the other hand the basic limitation of the "first 
Singular case'' (See Section V) is that in general the type 
number of the plant is not preserved, and for the "third 
Singular case" (See Section VII) is that the plant has to 
be at least type two. | 

For the above reasons, the "second singular case" was 
further considered and some extensions of its compensation 


structure were studied. 


Noe, 





B. SECOND SINGULAR CASE COMPENSATION STRUCTURES 


i The Basic Structure 


The basic compensation structure of the "second 


Singular case" is shown in figure (6-4) where: 


ares kee ed 7 


s‘(s+p,) (stp,)... (stp) 


represents the transfer function of the plant to be conm- 


pensated. 
H(s) = a (8-2) 


represents a gain adjustment, introduced as a negative 


feedback from the output to the input of the plant, and 


(S2zZ)) (StZ2. )etce(S+z. ) 
ee(s) = &Fo(s) = pends ee oi Oe, ee (8-3) 


s‘(s+p,) (stp,)...(stp,) 


represents the transfer function of a cascade compensator 
whose poles and zeros are identical with these of the 
plant and its Root locus gain is the second adjustment of 
the compensation structure. 

This compensation structure defined above, cor- 
responds to a singular system which was studied in Section 


VI as the "second singular case". 
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Its transfer function has the general form: 


m 
PA 
Ek | | (s+Z,) 
=) 


R(s) oN 2 = = = 4 
s | [sp eaks\| ltsep oT |estz,d+axl es+z.) 
k=l k=l fel 9) fen 





(8-4) 


The error constant Kx was evaluated and it 


-B 
was found to be: 


| | c. 


m 
_ 8 j=1 J 
K = = AR Sk ene (8-5) 
n 


sx-B 
I |p 


=] K 


aa 


The characteristic equation also of this basic 
Singular compensation structure was found to be in the 
general form: 


L+akF, (s)+8k 26s) = 0 (8-6) 


F 
§ 
From this basic structure six modified singular 


compensation structures were derived. They were analyzed 


as follows: 
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Z2. Modified Structures 


a. First Modified Structure 


The first modified singular compensation 
structure is shown in figure (8-1). It was derived from 
the previous basic structure and has been analyzed in 
Section VI. 

In figure (8-1) G(s) denotes the transfer 
function of the plant to be compensated and G..(s) the 
transfer function of the singular cascade compensator, 


which according to equation (6-18) its transfer function 


1S: 
m 
a[ | (stzj) 
Shee UD) —EE—EE (8-7) 
3N | |cs+p ) +a | | (S+z.) 
k=1 K j=l J 


where sN P, and z represents the poles at the origin, 


ale 
the other poles and the zeros of the plant respectively. 
The transfer function of this compensation 


structure was found to be: 








(8-8) 
m 
2 
TGs) wai rea 
S = = 
} R(s) ne n m m 
s | |cs+p.) tas | |¢sep ) | |cs+2.)+ek alee 
Kad ki 0S jar jal? 
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The error constant | was also found (See 


equation (6-85)) to be: 


K ae eee SK (8-9) 


where K denotes the error constant of the plant, before 
it was compensated by the singular cascade compensator. 
The characteristic equation of this compensa- 


tion structure will be in the general form: 


L+aFg (s) #8kF,“(s) = 0 (8-10) 


The basic advantages of this modified singular 
compensation structure with respect to the previous basic 


one, are the following two: 


(1) Simpler Structure. The compensation 
Structure is simpler because only a single cascade filter 
is required. Historically it is realized that series 


compensation has played a prominent role in system design. 


(2) Greater Error Constant. Comparison of 
equations (8-5) and (8-9), reveals that the error constant 
of this structure is greater than the error constant of 
the previous basic one by a factor equal to the Root locus 


gain k of the plant. 


Mors: 





b. Second Modified Structure 


The second modified singular compensation 
structure is shown in figure (8-2). It was derived from 
the previous one by changing the location of the singular 
compensator from the forward path at the feedback path. 
Of course this modified singular structure is no longer 
a unity feedback control system. 

For this compensation structure the transfer 


function was evaluated and it was found to be: 


is*T lee, THe, )+ak r | (s+2. )° 


T>(s)= £4S) re a a ae’ 


R(s) 
s [deena tos’ | lew AT foes )+8k a (sta, " 


(8-11) 


The characteristic equation also of this 


compensation singular structure was found to be: 
? 
1+G(s) G (s) = 0 
or 
ltaFy(s)+8k Fy“(s) = 0 (8-12) 


In order for steady state accuracy to be in- 
vestigated the following considerations were made: 


The error of the system is: 


E=R-C=(1-T,)R Cok 
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Introducing equation (8-11) into equation (8-13) 


yields: 


li i! Th mm 
ZN 
oT Newall lore pera [1 csr2,)° 
=] =] 


E(s)= < eH R(s) 


nN m mM 
2N 2 2 
sol sepa tas! | cso [| see, deen |] (see, 
eel hl K j=l J j=l J 





(8-14) 


When a type one plant was considered and a ramp 
input, i.e., N=l and RG Bviere and when the final value 
theorem was applied to equation (8-14), assuming a stable 
system, the steady state error was evaluated as: 
E(t)=B..(t) = 
tr 

For a physical interpretation of the above re- 
sult the following considerations were also made on the 
compensation structure which is shown in figure (8-2). 

For a type one plant, and a step input, at 


rd 


Q. RBRecause 


feeaay State # has €o6 be zero, 1.¢., Ess 


EF = R-C , it implies that R.,* Cog: Since Cae is a 
constant and the system is not a unity feedback system 
om #f C..- Therefore Oe £ R. 

When a type one plant was again considered 


Pad 


with a ramp input, at steady state E has to be constant, 
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“a 


Ae. 4 E,,~ Constant. Because E = R-C this implies that 
R = C5: Since Gee is a ramp and the system is not a 


unity feedback systen, oe # Gx Therefore C..* R, 
which implies that an infinite error EQ. must be expected. 
This result was also derived previously by analytical 
study of equation (8-14). 

Based on these remarks and results, the con- 
clusion is that such a System does not reproduce the input 
Signal at the output or in other words at steady state the 
output does not follow the input, which may be an advan- 
tage or a disadvantage depending on the specific applica- 
tion, and which can be eliminated as will be seen by the 
following modified singular compensation structure, 

When a type one plant was considered again and 
a step input applied, the steady state error was found to 


be: 


= (1- & 
geo) (ley) 


where A is the value of the step input, i.e., at steady 
state the output differs from the input by the value of 


ee 
Ss 


In order fcr all the above mentioned problems 
to be eliminated when this is desired, the following 
considerations were made on this modified singular 


Beructure. 
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When R was considered as the input signal which 
the output of the system was desired to reproduce, l1.e., 
to follow, then the error will be given by equation: 


E = R-C (8-15) 


By considering also an input which is not followed 
by the output to be Eee then according to the notations 
used: 


ON (8-16) 
Introducing equation (8-16) into equation (8-15) yields: 
E = R-T4R (SS 7); 


By assuming that the actual input to the system 
R is correlated to the desired input R by a function Y, 
ee : 


a 


R = YR (8-18) 
then equation (8-17) becomes: 


E = R-T, YR = (1-T,Y) R (8-19) 


Introducing equation (8-11) into equation (8-19) yields: 


Tern ) +(a-ky)s Tle, C Iesee, )+(Bk-oky) C Mese2,9 


E(s) & ee 


ee’ ) +as T lor, ) atc )+8k [Jory ) 


(8-20) 
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In order for the steady state error to be 
finite for a ramp input to a type one plant or zero for a 


step input, Y has to be defined as: 


16 


Introducing equation (8-21) into equation (8-18) yields: 


a 


R=. R (8-22) 


QlwW 


After the above result was obtained a third 


modified structure was derived and is analyzed as follows: 


c. Third Modified Structure 


The third modified singular compensation 
structure is shown in figure (8-3), where R is the input to 
the system which is to be followed at steady state by the 
output of the system. 

The transfer function of this modified struc- 


ture was found to be: 


n m m 
ae 2 
KN ls | |(s#25) tok eae) 
=] j=l j=} 


Nl Nn Mm Til 
2N 2 2 
so] |p tes’ | | spp [| (stz,)+0k 1 | ste.) 
k=1 k=1 j=1 J atl J 


Q|mD 


2 &s) _ 
T,(s)= ne 


(8-23) 
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The characteristic equation was also found to 
be in the general form: 


lta Fo(s)+8k Fe“(s) = 0 (8-24) 


By the same procedure as in the case of the 
Piest modified structure, the error constant Ki,-3 for 


this case was found to be: 


ee LS ae Ue (see) 





where Ky denotes the error constant of the plant before it 
was compensated. ‘For the examples which were worked out 
iHmosection Vi, the ratio oB/ {a -Bk} turns to be greater 
than unity. Therefore, when such is the case the error 
constant of the uncompensated plant will be improved when 
the plant is compensated. 

When this compensation structure was compared 
with the first modified structure it was found that they 
have the same characteristic equation but different error 
constant and transfer function. Therefore the stable 
region on the a-f8 plane will be the same in both cases. 
Also for the same operating point (a,8), although they 
will have the same characteristic roots, the residues of 
these roots will be different in each case due to the dif- 


ferent numerator of their transfer functions. Since the 
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dominancy of a characteristic root is a function of its 
residue, in addition to the distance of the location of 
this root on the s-plane from the imaginary axis, a dif- 
ferent degree of dominancy for the operating singular 


root must be expected in general in each case. 


For a specific plant and a specific operating 
point on the o-f8 plane a further investigation of the 
dominancy of the operating singular root can be performed 
either by a computer simulation of each compensation 
structure aril comparing their time response characteristics 
with those corresponding to the singular root second order 
model or by analytical evaluation of the residues cor- 


responding to the characteristic roots of each case. 


d. Fourth Modified Structure 


The fourth modified singular structure is. 
Shown in figure (8-4). It was derived from the first one 
by properly defining the feedback compensator in order for 
these two compensation structures to have the same open 


Hoop transfer function, 1.e.: 


G(s) 


ee (8-26) 
1+G(s) G..(s) 


G(s) Fee (Os = 
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from which the transfer function of the filter ca cs)) was 


evaluated to be: 


nl Tl m TM 
sN[ |csppis® | |eseppeal Its+2,-8 | |tst2,9) 
A.) = k=1 k=1 jr J jn 7 
sae a ; 
Bk [ ](s#25) 
j= 


(8-27) 


Because G(s) defined by equation (8-1) denotes the 
transfer function of a physical control system, it implies 
that N+n>m, i.e., the above defined filter has more zeros 
than poles. In practice a structure of this type is 
realized by introducing additional poles as needed, but 
placing them well outside the bandwidth of the system. 
Whether this can be done without affecting the singularity 
characteristics has not been investigated. Based on the 
argument requirement for singular conditions (See equation 
(3-12 and (10-11)), it is believed that if these additional 
poles are placed on the s-plane in such locations that 
their net contribution to the argument of G_.(s) is 
negligibly small, then it can be assumed that the singu- 


larity characteristics of the system remain unaffected. 


eal 





The transfer function of this structure was found 


to be: 


T,(s) = ey = G(s) 


ay : (8-28) 
1+ G(s)+G(s) G_ 


Gs) 


Cc 


or after some manipulations: 


Hitt 
2 
Bk | | (s+Z.) 


one ae a ae ea Ree aaa 
ON] [ese tasN | [espa [ I¢s+2,)+ ak[[es+2,9° 
k=1 k=1 j=1 j=1 
(8-29) 


The characteristic equation of this compensation 


structure was found to be in the general form: 


LtaFy(s) +8 k F,(s) = 0 (8-30) 


By the same procedure also, as in the case of the 


first modified structure, the error constant K_. 4 for 
this case was found to be: 

Ko =&8K (8-31) 
sx-4 a xX 


Comparison of the results which were obtained for 


the first modified compensation structure with these 


WAZ 





obtained above, reveals that the performance of the system 


in both cases is the same. 
e. Fifth Modified Structure 


The fifth modified singular compensation 
structure is shown in figure (8-5). It was derived from 
the second one shown in figure (8-2) by splitting the 
feedback path into two parallel paths as shown, such that 


their sum is the compensator transfer function: 
tt 
G..(s) = 1+G..(s) (8-32) 


from which the transfer function of the minor loop feed- 


back filter was found to be; 


k=1 
GG) = 


s. | | (stp, +0 | | (s+z,) 


k=1 j=l 


(8-33) 


m Nl | 
" : (o-adl, hesre,y-s | (S*+P,) 
Nn 


The transfer function of this compensation 


Seructure 1s: 


C(s) . @(s) 
14+G(s)+6(s) G(s) 








We 





or after some manipulations it becomes: 


qh m m 
rN] Hsp] |cst2j)t0x[ | cores)” 


T.(s) = tS} ee 


is 


Tony" +as T leer | ll bee )+8k | i(stz, )" 
j=l 


(8-34) 


The characteristic equation of this compensation 


structure was found to be in the general form: 
2 
l+aF,(s) 1 BkF,”(s) = 0 (8-35) 


Also when steady state accuracy considerations 
were made for this compensation structure the ores were 
the same as these for the second modified structure. Com- 
parison of the results corresponding to these two cases 
reveals that they are equivalent under performance con- 


siderations. 


f. Sixth Modified Structure 


This is shown in figure (8-6), and it was de- 
rived from the third modified structure with the same 
reasoning by which the fifth modified structure was de- 


rived from the second one. 
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It was also found that the third and sixth 
modified structures have the same transfer function, 
characteristic equation, error constant and thus under 


system performance considerations they are equivalent. 
3. Remarks 


Although the "second singular case" can be used 
for compensation of any control system, ‘Oo. it may have 
general application, its "basic structure" has certain 
disadvantages, which have been stated in the previous 
analysis and which may not be acceptable in some cases. 

When the characteristic equation of the basic 


structure is considered, it will have the general form: 
L+akF, (s)+8kF, “(s) = 0 


On the other hand the six modified structures, 
which were analyzed above, have certain advantages over 
the "basic structure’! and a common characteristic, i.e., 
they have the same characteristic equation, which has the 
general form: 

l+aF, (s)+8kF “(s) = 0 

g § 

Also their performance characteristics can be mathematical- 
inv ecorre lated . 

This fact aids the designer in making his decision 


on the problem of the structure without additional design 


effort and protects from arbitrary choices. 
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In the following example all these modified 
singular compensation structures were used to compensate 
the same plant in order that the derived analytical 
results could be compared. The compensated system per- 


formance was also studied in each of these cases. 


4. Example (8-1) 


In this example the plant of example (6-7) was 
considered for which the transfer function G(s) of the 


system before any singular filter has been introduced is: 


20(st2) 

G(s) = ———————_ (8-36) 

s(st1) (€st+5) (s+20) 

Also as operating point from the dominant root 
region of the a-8 plane corresponding to the singular root 
(c .=0.60, Was t-23); the point (a=920, §=2500) was ar- 
bitrarily selected. Then each of the six modified sin- 
gular compensation structures was considered separately 
and the performance of each for a step input, R=l1 Volt and 
a ramp input R=0.1 rad/sec, was studied in the computer as 


follows: 


a. First Modified Structure 


This is shown in figure (8-1). According 


to equations (8-7) and (8-36) the transfer function of 
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the singular cascade compensator G_.Cs) was found 


to be: 


2500 (st2) 
G_.(s) = 


By simulation of this system in the computer 
with the use of a DSL program the time response character- 
istics were obtained for a step and a ramp input. They 
are shown in figure (8-1A) and (8-1B) respectively. 


The characteristics of the transient response 


corresponding to the singular root (zc .=0.60, Wygrt-23) 
were evaluated and they were found to be: 

Time to peak of first overshoot: t = 3.19 sec 
Settling time oGey 5.41 sec 
Percent overshoot >: Mp = 9.48% 
Transient oscillatory frequency: w, = 0.99 rad/sec 


When the step input was considered then it 


was found from the computer print output that: 


Steady state output value : cme = R= 1 Volt 
Time to peak of first overshoot: th = 3.2 sec 
Settling time > t, = 3.4 sec 
Percent overshoot : Mp = 11% 
Transient oscillatory frequency: %, = 0.98 rad/sec 
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When the ramp input was considered then it 


was found also from the computer print output that: 


Steady state output value : Cae = Ro. 002 
Steady state output rate Co. = R= 0.1 
Settling time > ts = 5.4 sec 


These results can also be derived approximately from the 
corresponding figures of the time response characteristics, 


figure (8-1A,B). 


b. Second Modified Structure 


When the singular compensator G..Gs) was 
located at the feedback path as shown in figure (8-2), the 
transient response characteristics which were obtained are 
Shown in figure (8-2A) and (8-2B) for a step and a ramp 
input respectively. 

The exact transient and other characteristic 
values which were read from the computer print output have 


as follows: 


(1) Step Input. 


Steady state output value C= = 0.368 = = R 
Time to peak of first overshoot : ty = 3.2 sec 
Settling time Eg = Sot See 
Percent overshoot from one : Mp = 10.2% 
Transient oscillatory frequency ee ee) rad/ See 
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(2) Ramp Input. 


Steady state output value : Cae followed a ramp 


different than the input ramp. 


ae. 
ll 


Steady state output rate 0.036 = 


MR 


55 


5.4 sec 


Settling time cae 


c. Third Modified Structure 


When a gain adjustment equals with 2.717, which 


corresponds to the ratio §/a value, i.e: 


em 2.717 


AG, 





L[w 
1 


was introduced at the indicated in figure (8-3) location 
the third modified singular structure was obtained. 

The transient response characteristics which 
were obtained in this case for the same step and ramp 
input as in the previous cases, are shown in figure (8-3A) 
and (8-3B) respectively. 

The exact transient and other characteristic 
values which were read from the computer print output have 


as follows: 
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(1) Step Input. 


Steady state output value : Oe. = 1 = R 

Time to peak of first overshoot : t = 3.2 sec 
Settling time 6 ce Dt SCC 
Percent overshoot Mp = 10.2% 
Transient oscillatory frequency : %, = 0.92 rad/sec 

(2) Ramp Input, 

Steady state output value aren JOS OHO) 
Steady state output rate : Che = 0.1=R 
Settling time > ts = 5.4 sec 


d. Fourth Modified Structure 


The fourth modified structure which is shown in 
figure (8-4) requires a feedback filter Goacs) After the 
transfer function G_.(s) was evaluated from equations (8-27) 
and (8-36) the system was simulated in the computer. The 
transient response characteristics for a step R=1 Volt 
and a ramp R=0.1 rad/sec input, which were obtained in 
this case were the same with these corresponding to the 
first modified structure. They are shown in figure (8-4A) 


and (8-4B) respectively. 


180 





e, Fifth Modified Structure 


This modified structure involves a different 
feedback singular filter G_.(s) whose transfer function 
was evaluated from equations (8-33) and (8-36). The 
transient response characteristics for a step R=1 Volt 
and a ramp R=0.1 rad/sec input, which were obtained 
when this structure was simulated in the computer are 
Shown in figure (8-5A) and (8-5B) respectively, which are 
the same match these corresponding to the second modified 


structure. 


f. Sixth Modified Structure 


This is shown in figure (8-6) and was implemented 
from the previous one by introducing a gain adjustment as 
it 1s indicated in figure. 

Simulation of this system in the computer for 
a step R=1 Volt and a ramp R=0.1 rad/sec input, gave the 
Same transient and steady state characteristics as in the 
case corresponding to the third modified structure. They 
are also shown in figures (8-6A) and (8-6B) respectively 


for each input which was considered. 
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5. Comment 


The above presented results which were obtained 
from the computer not only confirm the corresponding 
analytical results but also make clear the flexibility 
which the designer has when compensation is performed by 


using singular line theory. 
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[X. DESIGN OF COMPENSATION BY USING A 


SINGULAR CASCADE COMPENSATOR 


A. INTRODUCTION 


Consider the plant which is shown in figure (4-1), 
and that this plant is going to be compensated by a sin- 
gular cascade compensator. The open loop transfer function 


of this plant is: 


(st+z,) (stz,)...(s+z)) 


G(s) = kF,(s) = k (9-1) 


s\(s+p,) (stp 2)...(stp,) 

It was assumed that the above transfer function is 
known and that the problem's specifications have been 
interpreted in terms of the desired dominant characteristic 
root (Zz, WW) for the compensated singular system. It was 
also assumed that tolerance on the polar coordinates ¢ and 
w, values for this root have been established from the 
problem's specifications. 

The compensated singular system is shown in figure 


(6-6) where G..(s) denotes the transfer function of the 


183 





required singular cascade compensator, whose standard form 
nS? 
(9-2) 


B(st+z_) (stz_)... (st 
ees): = 1 2 s+2) 
ae 


s\(stp,) (stp,)... (stp, )ta(stz,) (s+z,...(s+Z,) 


as was derived in Section VI C. 

Therefore the design problem in general will be limited 
to the problem of determining the required parameter pair 
values. For different types of applications this problem 
can be stated in terms of the parameter plane as: What must 
be the operating point, or the operating line segment, or 
the operating region of the parameter plane in order that 
the performance of the compensated singular system be as 
desired? 

The required design procedure in all these cases has 
been developed in Section VI. In this section the design 
process in the case where an operating point is desired to 
be defined on the parameter plane and from this the pre- 
cise form of the transfer function (9-2), was formulated 
in a sequence of basic design steps which must be followed 


mm a case. 


B. DESIGN PROCESS 


The required design process in general when an ac- 


ceptable operating point of the parameter plane is desired 
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to be found, is relatively simple and can usually be 
completed with satisfactory accuracy for most applications 
with graphical and mathematical techniques, by the fol- 


lowing four basic design steps. 


ie cLesSt otep 


The first step consists of drawing the Root locus 
of equation: 


kF,(s) = -1l (ges) 


2. Second Step 


The second step consists of insertion of an ini- 
tial cascade compensator, when the shape of the above Root 
locus is such that either an acceptable operating singular 
point cannot be located on this Root locus or although it 
can be located it is not acceptable under dominancy consid- 
erations (See VI, F,2). 

When an initial cascade compensator is required 
then this requirement will actually affect only the 
transfer function of the singular cascade compensator and 
neither the complexity of the rest of the design process 
nor the compensation structure (See example (6-6) and 
figure (6-28 )). 

In the simple case where such an initial compen- 
Sator will not be required then this step will consist of 
determining an acceptable operating singular point on the 


Root locus of equation (9-3). 
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3. Third Step 


The third step consists of determining the 
equation of the singular line corresponding to the 
selected operating singular point (See VI D3a equation 


(6-29)). 


4. Fourth Step 


This step which is the last one consists of de- 
termining the required operating point of the parameter 
plane, i.e., the required pair of a and 8 parameter values 


(See examples (6-6), (6-7), (6-8)). 
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ee EXTENSION OF THE SINGULAR LINE THEORY 


TQ PARAMETER SPACE 
A. INTRODUCTION 


The parameter space is an extension of the parameter 
plane when there are three or more parameters. 


Considering the characteristic equation: 
n 
F(s) = )- ay ah = 0 (10-1) 


where n 1S a real integer and ay 1s a real valued coef- 
ficient which may be a function of three or more parameters, 
i: @: 

a, = (es (ce Shee a eg nv!) (10-2) 


where net denotes the nth parameter. 


The point s can be defined in polar coordinates as: 


a jw, 4 i-c* GlO= 3) 


D. SINGULAR CONDITIONS 


If for a given pair ¢., w,, values, where |e. |<1 and 


ns 


oO equation (10-1) has an infinite number of real valued 
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Pee SEO 160. BS ¥'5-0 5, <2 .% nth parameter values which satisfy 


it, then for this equation singular conditions exist. In 
this case the point Cr We) is a singular point. This 
equation also for which singular conditions exist is 


called singular. 
ie. SINGULAR CONDITIONS ANALYSIS IN THE LINEAR CASE 


A simple case which was analyzed is the linear case, 
in which the variable parameters appear in the coefficients 
of the characteristic polynomial in the general form: 


ay. = by AFC B+dy oi OF ss F arene 


th 


+{Coefficient of the n parameters}, (n*? parameter }+ 


+{Constant term}, 
where Dy» Cho +++ are real constants. 
In this case equation (10-1) can be rearranged in the 
form: 
f, (s)atf£,(s) B+f£.(s) yt ...4£ (sin parameter}+ 
Pies) = 0 (10-4) 


where £, Gs), f(s), ... are polynomials of s. By analogy 
as in the case of the parameter plane (See Section IIB2) 


equation (10-4) can be written as; 


oJ nth, 
1 = 


je 59 j9 jo, th 
r,e 1 gtr,er °2 B+r ,€ 3 yt+...+re nin parameter }+r_ 0 
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or 


Bee ose sine; jatr, (cose, +jsine,)B+...+ 


Sean th wet: 
+r, (cose, +jsine,) {n parameter}+r_,, (cose, ,+jsine_,,) = 0 


Equating then real and imaginary parts to zero, 


yields: 


— 3 s th . 
Y,s1n@ oatr  sinée +...+ Y Sin n +7 sine = Q 
i ; ; 7 B = ea } nel nt) 


and 


il 
/ 


th 
+ ca See Ga 
r,Ccos@, a rcos8, B r,6086,,{ n~ } cosé 


n+} l 


According to [1] the above system has an infinite 


number of solutions if and only if: 


fe, 51n Moe r_si eeeor 
; 6,  T,81n6,_ Dan’ oy _ 7,51ne, _ The N64, 
i r 
r,cosé, r,cosé, Y,cose, noose, nto n+) 
or 
@, = = Hos n@. = tane 10-5 
tan6, = tang, = tane, tang, ae ( ) 


Equation (10-5) implies that for the characteristic 
equation (10-1) or (10-4) singular conditions exist if 


and only if: 
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> > 
Bh el 
— = 
ih ¢ 2) 
We a 7 
i} 
HA A 
=| =| 
+ + 
bs 4 
= = 
€ 2) ws 
Net Nel 
+ + 
A A 
Rh - 
= <= 





(10-6) 
/£2(s) £4108) + Kz O 
/£,, (8) £41 (5) elke 
Assuming that Eats) # 0 and dividing equation 
(10-4) by ere ee yields: 
fs) £2 Cs) f(s) f(s) 
ee ee. fy, neo! Se A ere 
£441 68) £i+1 0S) sigs sod aie ce 
(10-7) 
or 
F,(s) a+ F,(s) 6+ F,(s) y+...+ F.(s) n@@+1 = 0 
a 2 5 a n 
(10-8) 


When the characteristic equation is considered in the 
form of equation (10-8) then the necessary and sufficient 
conditions for singular conditions given by equation 


(10-6) become: 
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TJ 
— 
on 
77) 
ed 
il 
van 
— 
= 


(10-9) 


By 
le) 
rer, 
#2) 
a 2 
I 
wn 
ee) 
= 


where any of the K,, K ...K, is real integers. 


D. SINGULAR THEORY AND CONTROL SYSTEMS 


From equations (10-8) and (10-9) and by the same pro- 
cedure which was used in Section III for the case of the 
two parameters, it was found that in this case of the 
parameter space where there are n parameters in the 
characteristic equation of a control system the general 


form of a singular characteristic equation is: 


i m m m 
1 Z 
aC, F (s)+ 8C5F (s)+y C.F Dg emeae int }c F Hy eee =0 


1 

(10-10) 

where Ci; C,, Cay+-s; C> Chel are real constants and m4 
Mo, Mz,---, M, are n different real integers. (Otherwise 


when there are two or more of these identical then the 
number of the actual parameters can be reduced by combining 


terms of equal powers of F(s)). 
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For equation (10-10) the necessary and sufficient 


conditions, for singular conditions are given by equation: 


/F(s) = KI (10-11) 


where K is a real integer. 
Equation (10-11) implies that all the points of the 
Root locus of equation: 


KF(s) = -l Gilt) ay 


which do not tele on the real axis are singular points for 
the characteristic equation (10-10), where K in this case 
1s considered to be the Root locus variable parameter which 
takes values from zero to plus infinity and from zero to 
minus infinity. 

Finally any control system with n variable parameters 
appearing in its characteristic equation linearly depen- 
dent, is singular, if and only if, its characteristic 
equation can be written in the form of equation (10-10) 
which in the simplest case takes the form: 


wk 
n 


1+ oF(s) + @F°(s) + yF°(s) +...n°™ B@(5s) = 0 (10-13) 


1. Equation of a Singular n Dimensional Surface 


When an n dimensional parameter space was con- 
sidered then the general form of the corresponding singular 
characteristic equation of a control system is given by 


equation (10-10). Then according to the previous analysis 


| 





when the value of a singular point s. is substituted in 


equation (10-10) it becomes: 


ral 


i 


+ eC, + yC, t...t+n co +c = 0 (10-14) 


aC Z n el 


a A “ A “ 


where Ci> C,; Capes, e C+] are real constants. 


nN?’ 


Equation (10-14) is the equation of the singular n 
dimensional surface, corresponding to the singular point 


Ss which was considered. 


s? 
In order to determine a singular point of this 

Singular surface the values of n-1l parameters must be 

defined and then the value of the parameter left can be 

evaluated from equation (10-14). Thus one among an infinite 

number of real valued sets of parameter values can be de- 


fined, which in other words defines the coordinates of a 


Singular point in the parameter space. 


i. COMMENT 


From the previous analysis it can be seen that any 
non-linear case, either in the parameter plane on in the 
parameter space can be studied by defining any different 
parameter's product as a new variable parameter, e.g., the 
non-linear case in the parameter plane can be considered 


as a linear case in the parameter space (three dimensions 
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for this case) by defining the product a8 equals with y 
which will be the third variable parameter. 

After transforming a non-linear characteristic 
equation in a linear characteristic equation as it was 
described above then all these which have been stated for 
the linear case can be applied for the transformed 


characteristic equation. 
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XI. CONCLUSIONS AND RECOMMENDATIONS 


FOR FURTHER WORK 


This thesis has developed practical methods for de- 
Signing singular compensators, i.e., compensators which 
produce the desired Singular conditions in a single input- 
single output control system. 

The key to the development of these methods was the 
form of the characteristic equation which was considered, 
1 6fek 


af, (s)+8f£,(s) +f, (s) = 0 


where: f(s), f,(s) and f(s) are polynomials of s and 
a,8 are real variable parameters. 

Necessary and sufficient conditions for the above 
characteristic equation to be singular were determined and 
the singular line theory was then investigated to make 
possible the-development of the required standard set of 
rules for the design engineer. 

Application of singular line theory in control system 
design led to the conclusion that it can be used in 
general for design of compensation of any single input 
Single output linear control system with the following 


basic advantages. 


IESIS 





Required filters for any desired compensation 
structure have transfer functions which are known as a 
function of the two parameters a and 8g. The parameter 
plane diagram required for analysis and synthesis of a 
design problem consists of straight lines, i.e., the 
Singular lines which can be drawn by hand, without needing 
a computer. Design techniques and procedures are charac- 
terized by simplicity, flexibility and accuracy. Finally, 
tolerances on the parameter a and 8 values can be pre- 
cisely established. 

When considering the sensitivity of the transient 
response to parameter variations, it was found that the 
control system is less sensitive to variations in a and 8 
if the operating point on the singular line is remote from 
the point at which the singular line is tangent to the 
constant ¢, and Wig Curves. 

The basic design tools of the singular line theory 
design method are: The Root locus of the plant and the 
parameter plane diagram of the compensated singular 
System, 1.e., under certain conditions or in other words 
Singular conditions, a classical design method (The Root 
locus method), and a modern design method (the parameter 
plane method) can work together in such a way that each 


compensates disadvantages of the other. 
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For example stability analysis which is not a simple 
process in the parameter plane is performed on the Root 
locus diagram in the s-plane or studying the effect of 
variations in either one or both of the parameter a and 8B 
values on the transient response of the eres which is 
not a simple process on the Roots tocas diagram is per- 
formed in the parameter plane, etc. 

Conceptually the compatibility of the Root locus 
classical design method, which is adequate when there is 
a single variable parameter with the parameter plane 
modern design method, which is adequate when there are 
two variable parameters, when singular conditions exist, 
is based on the fact that in this case one of the para- 
meters can be defined through the equation of a singular 
line, as a function of the other. Thus a two parameter. 
problem becomes a one parameter problem. The similarity 
also of the Root locus of the plant with the Root locus 
of the compensated singular system, independently of which 
singular line is considered generalized the validity of 
the above concept and nade possible the development of 
graphical design techniques which simplify the work and 
markedly increase the power of visualization of the de- 
Signer. 

It was also realized that the singular line theory 
could be applicable in self-adaptive control system, since 


movement of the operating point (a,8) along a dominant 
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line segment of a singular line or within a dominant root 
region in the parameter plane implies almost fixed system 
performance. In such a case the adaptive controller need 
only to drive the system back to the dominant root line 
segment or dominant root region and not to a specific 
operating point in the parameter plane. 

The singular line theory was further extended in the 
parameter space and in the most general case where products 
of the variable parameters appear in the characteristic 
equation. 

Further study is needed for the case of a singular 
system with a negative error constant. This study can be 
initialized from an analysis of equation (8-25). 

Also the correlation, if any, of a singular feedback 
compensator with the linear state variable feedback maybe 
interesting. 

Some other related areas for future work, which are 
recommended are: Applications of the singular line theory 
in self-adaptive control systems. Singular surface theory 
and control systems, Singular line theory and multi-input, 
multi-output linear control systems, and finally any of 
the above in the general case where products of the variable 
parameters also appear in the characteristic equation, l.e., 


non-linear cases. 
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APPENDIX H 
SINGULAR POINT PROGRAM 


The purpose of this program is to determine the roots of a 
characteristic equation which are singular points with respect to the 
parameters qa and'8. 


This program is applicable to polynomials whose coefficients are 
of the form (b atc 8 +d), whereaandg are variable parameters and b, c, 
and d are real constants. 


For a given Characteristic equation this program solves the 
system (2-4), by the Cramer's rule method for a desired value of z 
and for a range of values for wn, which has been determined within 
the program and can be changed when it is desired, by changing the 
card 0035. 


The four columns which are generated at the print output are the 
following: The first is a list of the values for wn, which were used 
for the calculations and the other three are the corresponding nume- 
rators and common denominator from the Cramer's rule solution of the 
system (2-4). 


Any value of ®, in the first column, for which the associated 
values of the other three columns are zero's, together with the value 
of ¢ which was used defines a singular point for the characteristic 
equation. 


The input data for the program are: 
a. First card: 
Value of Gand N, format (F4.2, I2) where ¢ is the desired 
value of damping ratio and N is the order of the characteristic equation. 


b. Second card: 
Constant coefficients in ascending order (8E 10.5) 


c. Third card: 
Alfa coefficients in ascending order (8E 10.5) 


d. Fourth card: 
Beta coefficients in ascending order (8E 10.5) 
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APPENDIX I 


BEGURES: 





FIGURE (4-1) 
Block diagram of a unity feedback system 





FIGURE (4-2) 
First compensation structure 


UZ, 








FIGURE (4-3) 
Second compensation structure 





FIGURE (4-4) 
Third compensation structure 
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FIGURE (4-5) 
Fourth compensation structure. 
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FIGURE (4-6) 
Fifth compensation structure. 
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(5-1) 
Parameter plane diagram. 
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UNITS INCH 
UNITS INCH. 





X- SCALE = 2.00E+Ol 
Y~ SCALE 6.00E+r Ol 
Example (5-1) 
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X-SCALE =2.00 E+01 UNITS INCH. 
Y-SCALE =6.00 E+01 UNITS INCH. 
FIGURE (5-2) 
Example (5-1):Parameter plane diagram showing 
a singular line. 
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X-SCALE*5.00E-0! UNITS INCH. 
Y- SCALE =4 OOE+00 UNITS INCH. 


FIGURE (6-2) | 
Example (6-1):Parameter plane diagram. 
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FIGURE (6-3) 
Block diagram of the plant. 





FIGURE (6-4) 
Block diagram of the compensated singu- 
lar system corresponding to the “second 
Singular case",scheme. 
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FIGURE (6-5) 
Block diagram of the compensated singular 
system corresponding to the modified second 
Singular case,scheme. 
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FIGURE (6-6) 
Simplified block diagram of the modified 
second singular case,scheme. 
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FIGURE (6-7) 
Example (6-2): Root locus of the uncompen- 
sated system. 
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X-SCALE® 1LOOE+02 UNITS INCH. 
Y- SCALE= 1L.OOE+02 UNITS INCH. 


FIGURE (6-8) 
Example (6-2): Parameter plane diagram,K=20. 
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FIGURE (6-9) 
Example (6-2): Root locus of the compensated singu- 
lar system corresponding to the singular root: 
s=-0.3884+j0.8899 or (5,=0.40,%.,=0.97096) 
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X-SCALE 2 LOOE+O02 UNITS 
Y-SCALE * 1OOE+02 UNITS 
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FIGURE (6-10) 
Example (6-2): Parameter plane diagram,K=35. 
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FIGURE (6-11) 
(6-3): S-plane diagram. 
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FIGURE (6-19) 
Example(6-5): Root locus of the compensated singular 
system corresponding to the singular root (c¢,=0.30, 
Wye 2.15268). 
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FIGURE (6-24) 
(6-5): Root locus of the compensated sin- 


gular system corresponding to the singular root 
(c7,.=0.75,uns=0.60667). 
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FIGURE (6-26) 
Example (6-6): Block diagram of 
the plant. 





FIGURE (6-27) 
Example (6-6): Block diagram of the 
initially compensated plant. 
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FIGURE (6-29) 
Example (6-6): Root locus of the uncompensated system. 
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FIGURE (6-30) 
Example (6-6): Root locus of the plant after 
1t was compensated by an initial cascade con- 
pensated with transfer function(s+2)/(s+20). 
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FIGURE (6-31) 
Example (6-6):Root locus of the singular 
compensated system corresponding to the 
singular root (z =0 6,w, =I .2326) or (s= 
-0.7395+j0 .9861) 


249 





“yndut days e uojs warsds Ae_~nbuts payesuadwod 9a 


00 2r 00'hT 


SO ee PD ES a A a mm eincinAEGin unr uae sendy up or 


| 


00 'TE 00-8 


ss 
_ Ce a ee ee 


+ ee ee ee 


(Z€-9) J4YNdOI4 


00'S} 00°F 


ta 


Y2 $0 asSuodsayu 


(286°82=9 ‘G/=y -SONLPA AdJoOWeUeY) 
owti:(9-9) ayduexg 


"YOUT/SILUN OO T=a1e9S-4 
“YOUT/SItun OO°E=91e95-yx 


00°6 00°'9 


00°€ 


OOO} I- 


£C Om 
00°OF 


QO°O2e 


250 





"andut 


Ga°Let 


(9G8°6t=9 SE°TS=V >SONLPA AdJPoWeUed) 


dazs e@ uOof waashs vetnButrs pazyesuadwod ayy 4o asuodsau OWLI: (9-9) aLduex 


(€€-9) FYNOIS 


C0°hT 00'TE OD-8} 90 °S] 


"YDUT/SILUN OO'T=aLeoS-A 
"youT/Salun QO*E=aLeoS-x 


00 ‘TI 00°6 oD°'9S 00°& 


£C atxg 
t 


OO *Op 


aa *p- 


ao°O} 


70'o7 


0.0 "D& 


251 





(pyZ9°8/Z=4 SO6=V :Sanjea ADLIWeUL d) 
“yndut days e uoy wazysks uepnButs pazesuadwod 94} $O aSuodsau awWt{:(9-9) a,dwexy 


(y€-9) FyNOIS 


"YDUT/STLU OO*T=aLeos-A 
"YDUT/STLUN OO*E=31e9S-yX 


COLE Gorht oor OO'8r 00'S} Q0°Y) G0°6 00° 00's 


O27 Gi- 


— 

a Ss Se ey = pa caret 
——— ae me a ree er ee re a ey ne Sia Soe < P = a Bete 
es te ne ee en eee a ee ee ee ee re re eet a a 


| 


= C atx) 
00‘aF 


Q0°Or 


a0°Q& 


252 





(8S6/°TITl=9 ‘OOT=V :Sanjea Aa}Paweued) 
"yndut daqzs e& uot wazsks uetnBurs payesuadwods auy 40 asuodsau SWLI?(9-9) aLduexg 


(S€-9) 3YyNOI4 
"YOUT/SILUN OO T=8Le9S-A 
"YOUT/S2LUN OO*E=aLe9S-YX 


63 Le 00°ht Oa OD 'B} 00 ‘SF 00 ‘GF 00°t 00°9 00°€ 


i “eae a a a re a a a A A ee a 


* 


C ox) 


t- 


Q0°OP- 


00°OFr Q'0°Qt 


00 °@E 


253 





(ZvO'SZT=9 SOZI=V :Sanzea AapawWeued) 
"yndut days e uot waqsks ue_nButs payesuadwos auzy Jo asuodsau awit: (9-9) atdwex 


(9€-9) JYNOIS 


"YOUT/STLUN OO T=aLe9S-A 
"YOUT/STLUN OO E=aLe9S-xX 


0O*LT G0°h?’ 00°Te GO° BI OD°S} 00°C} OO0'6 OD’d 00°E 
Me 
9 
o 
oS 
vee Le) 
of 
Zh 
gf | od 
a a ae eee CE : “ 9 
{ oe 7 = 3 
ad 
L 


elope? a 


20°V@E 





00°4LU 


(Ilvpl‘'9IZ=G SG°IELTHaY :San[ea AdzTowWesed) 
"gyndut dazs Be uoy waqgsks ue[tnbuts payesuadwod ay. Jo asuodsau awi{:(9-9) a,dwex 


(/€-9) JYNOI4 


"YDUT/STLUN OO T=aLe9S-4 
"YOUT/stlun OO" E=91e9S-¥X 


OO'he 00°Te DD SF 00 ‘ST OD°eT DO °6 po "9 OD°E 


a OR, 4 
- 


{[ OX? 


Tt 


Oo ‘O- 


00’02 Q0°‘OF 


OO * Ce 


Ess 











-andut 


OLE 


A 


(prO'€pp=d SOOZ=y :Sanpea AazTaweued) 


dazs e wos Waysks ue[nBbuts payesuadwos auz 40 asuodsau 


(Sieso)) ation 


ap hy OO°Te 0 ‘Bt 00°Sr G0°Ur 


SS TS A I OS SS TT ED aS SS ET NE UE A ee Neh yl a we wn ee 


owL{:(9-9) a ,dwex 


"YDUT/SFLUN OO" T=aLe9S-A 
"YOUT/SFLUN OO" E=9Le9S-X 


00°G 00°9 


ee 
SS Re a a Do Re re 


00°E 


O0G°ON= 


@o'ag 


Q0°CE 





(Z99°E€/8=0 SOEEHV :SANLeA AaQoweUueg) 
"andut days e uoj waqgsks ueypnbuts payesuadwod ayy jo asuodsau OWwLi: (9-9) aLdwexy 


(6€-9) JYNOT4 


00-°LT OO°HE OO'Te Oo-gF 00-S] 00 ‘TF 


nn ee ee te rt as nec i mn 


"YOUT/S2LUN OO'T=ALe2S-A 
"YOUT/S2LUN OO E=aLPIS-¥ 


00'6 da:9 00°€ 


Qa ‘Oy 


QQ°Oe Qd°pP 


ODOR 





(VES°SOTT=9 SOObP=V :SANLPA UaQZoWeUued) 
“yndut days e uwoy WazskS uetNbuLs paresuadwod ayz JO asuodsau awiy:(9-9) a, dwexy 


(Ov-9) JYNdIS 
"YOUT/SILUA OO’ T=aLeoS-A 
"yYoUT/STLUN OO? E=9Le9S-¥X 


G0°40 OO 'hT O0'Te 00 “BI 00 °S} OO‘EF Q0 ~6 00°5 00°¢ 


A teen nes sre = 
SE PS SS CT LS AT SO SS SS SS A oS a a A NR SE A | A a GPS SS el Pe og ne gee 2 een t~ a ee, 
{ Sarenennen 


20 oO fo 


ao ‘a? Qa OF 


DO" 


298 





(012°660Z2=9 SOOL=V :SanleA UazZoaweued) 
"ayndut daqzs e@ uoj wazsAs uetnbuts payesuadwod ayy Jo asuodsau auty]:(9-9) aL duex] 


(Ip-9) JYNdIS 
“YOUT/SILUM OO’ T=3LP2S-A 
“YoUuT/STLUN OO*E=3Le9S-X 


OODLE GO-Phe 6O'1E S08} OO°Ssy 00 SF 00°6 G0°9 ao 'f 


~~. 


——— a 
{ na ee re eS a EY I I I A a I 2-2, 


a eat 


Ye 


, 


OG°D) 


GQ°O— 


£ atx 
ole ir) 2 


T- 


Qeoe 


00 °0E 


250 








(900° €60€=9 SOO0T=V :san,teA uazaweueg) 
"yndut days e voy WaSks uepnBurs payesuadwod 94} JO aSuodsad awit: (9-9) aptduexy 


(Zv-9) JyUNdI4 


"YOUT/SILUN OO T=aLe9S-4 
"YOUT/SILUN OO'E=aLe9S-xX 


00L”e OohT Oo 'T? 00 °8} 90°51 oleh | go'6 09°9 O0°E 


OS A SO EE EL ES SS I ST NT SS A i A EE A aS re Ny nh SP a eh pay A he ah er eS es a ae A A A RS SL I OY a “ 


—_ 


[ ax] 


@Q> Qu- 


G0°OF 


Qa° og 


0O°0e 


260 





“yndut days e uoy wazysks uepnBults paryesuad 
(€7-9) 


Q0°LE OO°HE’ 0g ‘tt O0-§} 00°57 


IP A Sl a ee —_—S eee —_—_eooe -- — ~ ———_—————— —_ —— Sl i PD a a ie OD: 
_ ee > ew Sos ee —_ — See eae ee a - — 


(pL°980v=4 SOOET=¥ :santea uaqyaueueg) 
WOD 94} JO asuodsau aL: (9-9)a;duex] 
JUNSI4 


"YOUT/SILUN OO T=aLeos-4 
"YDUT/SILU OO E=aLeoS-X 


OO'EF OD'é 00°9 O0°E 


Qp°Dd _: 


ole bia) ag 


ao*at 


261 





(€2°6tLb=9 SOOSGT=V :san[ea sazaweueg) 
"gndut days e vos waqsks uetnburs pazesuadwods ay jo asuodsau oWt{:(9-9) a ,dwex 


(vv-9) 3YNOI4 
“YOUL/STLUQ OO’ T=aLe2S-A 
"yoUT/STLUN OO’ E=3aLe9S-X 


C0'Lz O0°Ae C0°st 00 °8? OO'SF OO Kw! 00°6 00°93 O2°'¢e 


0@°O}r- 


°“O- 


Ga uy” 


O0O°0E 


262 





(SG8°vI6-=9 SOSST=V :Ssanyea Aazoaweueq) 
"3ndut days e uoj waqsks uernBurs pazyesuaduon 942 $0 asuodsad awtt:(9-9) atdwexy 


(Sv-9) JUuNOIS 


"YOUT/SPLUN OO*T=aLeoS-A 
"YoUT/SztLUN OO E=aLeos-y 


~GO *O- 


00° C0°he OO'N? OD St Oo-S} eo°tT 00'S g0°9 00° 

‘ 
S 
| 
o 
ue mr re 
PO a aa fal 
ot ee a 3) 

tel 

Am 


20‘at 


G.0°ae 


263 





(22L°IlvS=d SOOZT=V :San[eA uazoueued) 
"yndut dazs e wos waqsks ueynBbulrs paqesuadwod auz 40 asuodsau 9WLI:(9-9) aL dwex] 


(97-9) JYNOIS 


"YOUT/STLUN OO*T=8L289S-) 
"YOUT/S2LUN 00*E=a1e9S-¥ 


QODLE OHE O0°Fe CO'sr O0°S} 00'°GF Q0 °6 00°9 00°E 


~~ a 
a SS AY SAP A ESSE SY SEEN RO, tee SY ers UNE aL GE es SE-B. mS OEeE eh eeneEeEd aagv ome ee a ert ET BS a ty ge TC a ae tt i en OO ae as “> a J 
- 


eee ees 


c os 


eee 


04° Op- 


OO°OF 


QQ" OF 


264 





(896°2pLG=4 ‘OO8T=¥ :Sanlea sazameued) 
-yndut days e uoy WaysAs Uue_NBults payesuadwod ayy JO asuodsau ality: (9-9) aL dwexy 


(4-9) JYNOIS 


"YOUT/STLUN OO*T=9Le9S-h 
"YOUT/STLUN OO" E=3LP9S-X 


cole OO'hE COT OO * 3} Oo°sT OO'UF op°G o0o'9 a0°e 


20°O]- 


— 
en A RY A ar a gy A A ger el gy, are a ag eR pet yt a ge AP ae S, ie 
—— 


00°0e 
{ 


Q0°0E 








(Sth°SOv9=d *O00Z=V :Ssanjea AVTPOWeUR  ) 
"gndut daqzs e voy waqsks uetnburs payesuadwoo 942 $0 asuodsau awit: (9-9) a, duwexy 


(8~-9) JYNOI4 


"YOUT/SZTLUN OO" T=aLe9S-A 
"YOUT/SZLUN OO*E=aLe9S-yX 


00e Ar gore OO'st 00's} 00°87 03°65 00'9 OO'Ee 


_— 


Nee 


— 


LD VIT LVAD ADP DYPVAP DAA PDI ISVS I 


5- Carn] 


99 


O0°OF 


Oo" O1/—- 


“Oo 


@o-ae 


Go°0k 


266 








o/0 


‘ 
& 
Pe 
S 
3 * 
% 
tf 
My 
4 
f 
(° 
\o wate \ 
: wre 
wor 
if » 
oy ye 
6 ; 
S = 
) Sas 
Ce, 
fe, 
2, 
> 
=, 
° 
By 
S 
o 
9 
3 
000 005 olo Os o2O ox 
3 
° 
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FIGURE (6-49) 
Example (6-6): Parameter plane diagram of the 
compensated singular system.Singular point 
(7, =0.6 07 5=1.232586) 
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FIGURE (6-50) 
Example (6-6): Parameter plane diagram of the 
compensated singular system.Singular point 
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FIGURE (6-51) 
Example (6-7):Root locus of the compensa - 
ted singular system. 
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FIGURE (6-52) 
Example(6-7):Parameter plane diagram and 
boundary lines. 
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FIGURE (7-1) 
Example (7-1): Block diagram of the 
initial system (Plant). 
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FIGURE (7-2) 
Example (7-2): Block diagram of the 
compensated singular system. 
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FIGURE (8-1) 
First modified singular compensation 
Structure. 





FIGURE (8-2) 
Second modified singular compensation 
Structure. 
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FIGURE (8-3) 
Third modified singular compensation 
Structure. 





FIGURE (8-4) 
Fourth modified singular compensation 
structure. 
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FIGURE (8-5) 
Fifth modified singular compensation 
) structure. 





FIGURE (8-6) 
Sixth modified singular compensation 
Structure. 
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